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We analyze pairing of fermions in two dimensions for fully-gapped cases with broken parity (P) 
and time-reversal (T), especially cases in which the gap function is an orbital angular momentum 
{I) eigenstate, in particular I = —1 (p-wave, spinless or spin-triplet) and I — —2 (d-wave, spin- 
singlet). For I 7^ 0, these fall into two phases, weak and strong pairing, which may be distinguished 
topologically. In the cases with conserved spin, we derive explicitly the Hall conductivity for spin 
as the corresponding topological invariant. For the spinless p-wave case, the weak-pairing phase 
has a pair wavefunction that is asympototically the same as that in the Moore-Read (PfafBan) 
quantum Hall state, and we argue that its other properties (edge states, quasihole and toroidal 
ground states) are also the same, indicating that nonabelian statistics is a generic property of such 
a paired phase. The strong-pairing phase is an abelian state, and the transition between the two 
phases involves a bulk Majorana fermion, the mass of which changes sign at the transition. For 
the d-wave case, we argue that the Haldane-Rezayi state is not the generic behavior of a phase but 
describes the asymptotics at the critical point between weak and strong pairing, and has gapless 
fermion excitations in the bulk. In this case the weak-pairing phase is an abelian phase which has 
been considered previously. In the p-wave case with an unbroken U(l) symmetry, which can be 
applied to the double layer quantum Hall problem, the weak-pairing phase has the properties of the 
331 state, and with nonzero tunneling there is a transition to the Moore-Read phase. The effects of 
disorder on noninteracting quasiparticles are considered. The gapped phases survive, but there is 
an intermediate thermally-conducting phase in the spinless p-wave case, in which the quasiparticles 
are extended. 



I. INTRODUCTION 

Most theories of superconductivity, or more generally 
of superfluidity in fermion systems, depend on the con- 
cept of a paired ground state introduced by Bardeen, 
Cooper and Schrieffer (BCS) in 1957 |J|. The ground 
state may be thought of loosely as a Bose condensate 
of pairs of particles, since such a pair can be viewed as 
a boson. Within BCS mean field theory, such a state 
forms whenever the interaction between the particles is 
attractive. For weak attractive interaction the elemen- 
tary excitations are fermions (BCS quasiparticles), which 
can be created by adding or removing particles from the 
system, or in even numbers by breaking the pairs in the 
ground state, and the minimum excitation energy oc- 
curs at fermion wavevector near kp, the Fermi surface 
that would exist in the normal Fermi liquid state at the 
same density of particles. There is also a collective mode, 
which is a gapless phonon-like mode in the absence of 
long-range interactions between the particles. This mode 
would also be present if one considered the pairs as ele- 
mentary bosons, and would be the only elementary low- 
energy excitation in that case. If the attractive interac- 
tion becomes strong, the energy to break a pair becomes 
large, and at all lower energies the system behaves like a 
Bose fluid of pairs. In the original BCS treatment, each 
pair of particles was in a relative s-wave (I = 0) state, and 



the minimum energy for a fermion excitation is then al- 
ways nonzero. No phase transition occurs as the coupling 
strength is increased to reach the Bose fluid regime. 

Not long after BCS, the theory was generalized to 
nonzero relative angular momentum (l) pairing, and af- 
ter a long search, p-wave pairing was observed in He"^ 
B. It is believed that d-wave pairing occurs in heavy 
fermion and high Tc superconductors. Some nonzero I 
paired states generally have vanishing energy gap at some 
points on the Fermi surface (for weak coupling), while 
others do not. While the absence of a transition is well- 
known in the s-wave case, it seems to be less well known 
that in some of these other cases, there is a phase tran- 
sition as the coupling becomes more strongly attractive. 
One reason for this is that the strong-coupling regime 
must have a gap for all BCS quasiparticle excitations. 
But even when the weak coupling regime is fully gapped, 
there may be a transition, and these will be considered 
in this paper, in two dimensions. 

In the paired states with non-zero / there are many 
exotic pheneomena, especially in the p-wave case, due to 
the breaking of spin-rotation and spatial-rotation sym- 
metries. These include textures in the order parameters 
for the pairing, such as domain walls, and quasiparticle 
excitations of vanishing excitation energy on these tex- 
tures (zero modes) (these are reviewed in Ref. ^J). In 
transport, there may be Hall- type conductivities for con- 
served quantities, such as spin and energy, which are pos- 



sible because of the breaking of both parity (P) and time 
reversal (T) symmetries. The breaking of these symme- 
tries, and topological aspects of the paired state, are more 
crucial for the ocurrence of these effects than is the angu- 
lar momentum of the pairing; the pairing need not be in a 
definite angular momentum state. Many of these effects 
have been discussed in remarkable papers by Volovik, of 
which a few are Refs. |§-||]. These are related to effects 
we will explore in this paper. We will make an effort to 
separate the effects related to breaking continuous sym- 
metries spontaneously, which leads to familiar Goldstone 
mode physics, from those connected with topological ef- 
fects, quasiparticle zero modes and Hall-type responses 
for unbroken symmetries. 

In this paper we will make extensive use of the meth- 
ods for BCS paired states, and consider the transitions 
between the weak and strong coupling regimes in two 
dimensions. In the weak-coupling regime, exotic phe- 
nomena are possible when parity and time reversal are 
broken. The results are applied to the fractional quan- 
tum Hall effect (FQHE) by using the composite fermion 
approach, to be reviewed below. We also consider effects 
of disorder on the phases and transitions, also within 
BCS mean field theory. In each Section, we try to make 
the initial discussion general and accessible for workers 
in many-body theory and superconductivity, before spe- 
cializing to applications to the FQHE. In the remainder 
of this Introduction we will give an overview of the back- 
ground and of the results of this paper. 

We now review some background in the FQHE P|. 
The original Laughlin states pffl occur at filling fac- 
tors V = l/g, with q odd (the filling factor is defined 
as i^ = n^o/B, where n is the density in two dimensions, 
B is the magnitude of the magnetic field, and $o = hc/e 
is the flux quantum) . An early idea of Halperin |ll|| was 
to generalize the Laughlin states by assuming that under 
some conditions, the electrons could form pairs, which as 
charge-2 bosons could form a Laughlin state. A variety 
of such states were proposed pi|-|l4| . Since the Laughlin 
states for bosons occur at filling factors for the lowest 
Landau level (LLL) i^b = 1/"^, with m even, and the 
filling factor for the electrons is related to that for the 
bosons by t^ = Aiyt, the electron filling factor is either of 
the form v — 1/q {q an integer), or v = 2/q, q odd. In 
particular, these values include f — 1/2, 1/4, . . . , which 
do not correspond to incompressible states in the usual 
hierarchy scheme [Q, because the filling factors in the 
latter always have odd denominators (when common fac- 
tors have been removed). 

The relation of the paired states in the FQHE to those 
in superfluidity theory becomes much closer once one in- 
troduces the notion of composite particles flfl-pfl]. A 
simple, direct formal approach is to use a flux attach- 
ment or Chern-Simons transformation (see e.g. Ref. pi] 
in particular), which represents each particle as (in the 
case of most interest) a fermion plus an integer number 



4> of (5-function flux tubes. After the transformation, the 
system can be represented by an action that includes a 
Chern-Simons term for a U(l) gauge field, that couples 
to the Fermi field. We refer to this as the CS fermion 
approach. The net magnetic field seen by the fermions 
is the sum of the external field and the (5-function fluxes 
on the other particles. In a mean field treatment, given 
a uniform density of particles, this produces a net aver- 
age field that vanishes when the filling factor is i^ = l/cj). 
In this case the fermions can form a Fermi sea pq |, or 
they could form a BCS paired state. Some of the exist- 
ing paired FQHE trial wavefunctions can be interpreted 
this way, as pointed out in Ref. |1J], and others can be 
constructed by analogy flj]. 

A more physical way of looking at the formation of the 
composite particles, particularly when they are consid- 
ered as the elementary excitations of the system, is as 
bound states of one of the underlying particles (or parti- 
cles for short), and 4> vortices in the particle wavefunction 
pi| , p6| . The bound states, which correspond to the CS 
fermions, again behave as fermions in zero net field if the 
particles obey Fermi statistics and is even, (or if the 
particles obey Bose statistics and cj) is odd) and v = !/(/). 
Note that we will consistently use the term "particle" 
for the underlying particles, and "fermion" for the CS 
or composite fermions (bound states). Some statements 
apply also when the transformed particles are bosons (ob- 
tained by interchanging the words "even" and "odd" in 
the preceding definitions), in which case we refer to com- 
posite particles. It is generally more important to keep 
track of the statistics and net magnetic field seen by the 
composite particles than those of the underlying parti- 
cles. Recent work has formalized the bound state pic- 
ture, and improved our understanding p7|-p0[. However, 
the results of the CS approach remain valid, and because 
that approach is simple to use, and we will mainly require 
only a mean field picture here, we assume that that is the 
approach we are implicitly using. 

The Laughlin states can be viewed as Bose condensates 
of composite bosons in zero net magnetic field ||l^,^,^ . 
Because the bosons are coupled to a gauge field (in the 
CS approach, the CS gauge field), vortex excitations cost 
only a finite energy, but there is still an effective Meissner 
effect for the CS gauge field. Because the fiux of the CS 
gauge field is related to the particle density, a vortex 
carries a fractional charge and corresponds to Laughlin's 
fractionally-charged quasiparticles jl^ (we refer to such 
excitations as FQHE quasiparticles). Hence the Meissner 
effect in the superfiuid becomes the incompressiblity of 
the FQHE state (there is of course no Meissner effect or 
superfluidity in the response to the electromagnetic fleld). 
Similarly, when pairing of composite fermions occurs in 
zero net magnetic field, the state becomes incompressible 
|lj|. In contrast, the Fermi liquid state of Halperin, Lee 
and one of the authors pq] has no Meissner effect for the 
CS field and is compressible pSl. 



The wavefunction proposed by Haldane and Rezayi 
(HR) in Ref. |]13[ is a spin-singlet paired state, which 
can be interpreted as a d-wave paired state of compos- 
ite fermions ||lj]. Moore and Read (MR) jlj] proposed a 
p-wave paired state (the "Pfafhan state" ) of spinless elec- 
trons with a structure analogous to the HR state. Both 
states can occur for filling factors v — Xjl^ 1/4, .... The 
HR state was proposed as an explanation for the observed 
V = 5/2 QH state |Q, which collapses when a parallel 
component of the magnetic field is applied, suggesting 
that it is a spin singlet. However, it was also proposed 
later that the 5/2 state is the MR state p^. In both 
proposals, the LLL is filled with electrons of both spins, 
and the paired FQHE state describes only the electrons 
in the first excited LL. The latter proposal is supported 
by recent numerical work [^,Q. The collapse of the 
state under a parallel magnetic field must then be due 
to another mechanism, involving the effects of the finite 
thickness of the single-particle wavefunctions in the direc- 
tion perpendicular to the two-dimensional layer, which is 
poorly understood at present. Another paired state with 
a similar interpretation is the 331 state |ll[ , which can be 
viewed as a p-wave paired state of two-component com- 
posite fermions ||l^,^,Q. It is likely that this is closely 
related to a FQHE state observed in double-layer and 
single-thick-layer systems at v = 1/2 [ p7| , p8[ . 

Moore and Read llj] suggested that nonabelian statis- 
tics might occur in QH states, and the Pfaffian state was 
proposed as an example. Nonabelian statistics means 
that the space of states for a collection of quasiparticles 
at fixed positions and quantum numbers is degenerate, 
and when quasiparticles are exchanged adiabatically (for 
which we need the system to have an energy gap for all 
excitations), the effect is a matrix operation on this space 
of degenerate states. This generalizes the idea of frac- 
tional statistics, in which the effect of an exchange is a 
phase factor and the states as specified are nondegener- 
ate; when the phase for an elementary exchange is ±1, 
one has bosons or fermions. The arguments in MR that 
this would occur were based heavily on the identification 
of the many-particle wavefunctions in the FQHE as chiral 
correlators (conformal blocks) in two-dimensional confor- 
mal field theory, which possess similar properties under 
monodromy (analytic continuation in their arguments). 
It was expected that an effective field theory description 
of these effects would be based on nonabelian Chern- 
Simons theories, which are known to be connected with 
conformal field theory, and to lead to nonabelian statis- 
tics |B3]. In the MR state, and other paired states, there 
are, apart from the usual Laughlin quasiparticles which 
contain a single fiux quantum, also finite energy vortex 
excitations containing a half fiux quantum p4| , and it is 
these which, in the MR state, are supposed to possess 
nonabelian statistics properties. 

Evidence for nonabelian statistics in the MR (Pfaf- 
fian) state accumulated in later work EQ-H, which in- 



vestigated the spectrum of edge states, quasihole states, 
and ground states on the torus (periodic boundary con- 
ditions), all of which were obtained as the zero-energy 
states for the three-body Hamiltonian of Greiter et at, 
for which the MR state is the exact unique ground state 
p2[ . The states found agreed precisely with the expec- 
tations based on conformal field theory. There was also 
evidence for similar effects in the HR state M,E3,H, 
however the interpretation was problematic because the 
natural conformal field theory for the bulk wavefunctions 
is nonunitary and therefore cannot directly describe the 
edge excitations, as it does in other cases such as the 
MR state. Some solutions to this were proposed |{42| , [l5| . 
Explicit derivations of nonabelian statistics and of effec- 
tive theories have been obtained later for the MR state 
[|6|,|3. The 331 state ^ is an abehan state, which 
can be viewed as a generalized hierarchy state |48| , [49| , 
as is evident from the plasma form of the state, and 
these two descriptions are related by a bosonization map- 
ping |42 4J|. The hierarchy states and their generaliza- 
tions possess abelian statistics properties, which can be 
characterized by a (Bravais) lattice |^,^ . Thus incom- 
pressible FQHE states in general can be divided into two 
classes, termed abelian and nonabelian. It is clear that in 
Halperin's picture |ll[ of bound electron pa.iis which form 
a Laughlin state of charge 2 bosons, the properties will 
be abelian and are simply described by a one-dimensional 
lattice, in the language of Ref. Esj (they are the simplest 
examples of a class of abelian states in which the ob- 
jects that Bose condense contain more than one electron 
plus some vortices, while the hierarchy states, and all 
generalizations considered in Ref. Esf, have condensates 
involving single electrons). 

In spite of the work that has been done, one may still 
ask questions such as what is the microscopic mecha- 
nism, in terms of composite fermions, for the degenera- 
cies of FQHE quasiparticle states that is the basis for 
nonabelian statistics, and whether it is robust against 
changes in the Hamiltonian. A similar question is about 
the effects of disorder. These will be addressed in this pa- 
per, by a direct and simple analysis of the paired states 
using BCS mean field theory, and developments such 
as the Bogoliubov-de Gennes equations ^^. We find 
that the nontrivial paired FQHE states are related to 
the weak-coupling regime (or more accurately, the weak- 
pairing phase); in particular the MR and 331 states have 
wavefunctions that contain the generic long-distance be- 
havior in spinless and spin-triplet p-wave weak-pairing 
phases, respectively. In contrast, the strong-coupling 
regimes, or strong-pairing phases, lead in the FQHE 
to the Halperin type behavior. There is a similar pic- 
ture in the spin-singlet d-wave case, except that the 
HR state, which might have been expected to represent 
the weak-pairing phase, is in fact at the phase transi- 
tion, and therefore is gapless in the bulk. The weak- 
pairing phases are topologically nontrivial, and possess 



edge states and nontrivial quantum numbers on the vor- 
tices (FQHE quasiparticles) ; the spinless p-wave case 
(MR phase) is nonabehan, while the spinful p-wave case 
(with unbroken U(l) symmetry), and the spin-singlet d- 
wave case are abelian states, which we characterise. We 
also consider the effect of tunneling on the double layer 
system which is the best candidate for realising the 331 
state, and show that the phase diagram includes a MR 
phase. The theory also leads to a description of the crit- 
ical theories at the transitions, at least within a mean 
field picture. The role of fluctuations, and the full ef- 
fective field theories at these transitions, remain to be 
understood. 

We also consider disorder within the same approxima- 
tion, making use of recent results on noninteracting BCS 
quasiparticles with disorder [pl|-p8| , and in particular we 
find that in the spinless p-wave case, there can be an in- 
termediate phase with the thermal properties of a metal, 
between the two localized phases which correspond to 
the weak (MR) and strong-pairing phases of the pure 
case. A disordered version of the MR phase still exists in 
the presence of disorder, though its properties, including 
nonabelian statistics, may become more subtle. 

One further issue which we discuss is the transport co- 
efficients of Hall type for various conserved quantities, es- 
pecially spin and energy. Concentrating on the quantities 
that are related to unbroken symmetries, we derive ex- 
plicitly the values of these conductivities in the spin case, 
for spin-singlet and triplet states, and show that they are 
quantized in the sense of being given by topological in- 
variants (in the disordered cases, we do not prove this 
directly). There have been claims that, in some sense or 
other, the ordinary Hall conductivity for charge (particle 
number) transport takes a nonzero quantized value in a 
He^^ film in the A phase ^^ and in a dx2 y2+idxv (i-e. P 
and T violating d-wave) superconductor [ p9y60| . It seems 
unlikely to us that these claims are correct, if the Hall 
conductivity is defined in the usual way, as the current 
response to an external (or better, to the total) electric 
field, taking the wavevector to zero before the frequency. 
While one can set up a detailed calculation, using a con- 
serving approximation as in the Appendix, which duly 
includes the collective mode effect in this case where the 
symmetry corresponding to the transported quantity is 
broken, we prefer to give here a more direct and appeal- 
ing argument. This works in the case where pairing is 
assumed to occur in a system of interacting fermions of 
mass m with Galilean invariance in zero magnetic field, 
as in most models of pairing. If one considers the linear 
response to an imposed uniform but finite frequency (lo) 
electric field, then there are well-known arguments that 
the conductivity is simply a^^{uj) = ne'^Sfj_^/{muj+ir]) {^, 
u = X, y, and 77 is a positive infinitesimal). This arises 
from the so-called diamagnetic term, and the contribu- 
tion of the two-point current-current function vanishes in 
this limit. This is independent of interactions, and hence 



also of whether the interactions produce pairing or not. 
The result can be understood as the contribution of the 
center of mass, which is accelerated by the applied uni- 
form electric field, while the relative motion of the parti- 
cles is unaffected, as a consequence of Galilean invariance. 
There is clearly no Hall conductivity. (However, the sim- 
ilar calculation in a magnetic field produces the standard 
Hall conductivity, and is an aspect of Kohn's theorem.) 
We may be curious about non-Galilean invariant mod- 
els, and whether a paired state could have a quantized 
Hall conductivity as claimed. But if it were quantized, 
it would be invariant under any continuous change in 
the Hamiltonian that preserves the gap. Hence, in the 
ground state of any model that can be continuously con- 
nected to the Galilean-invariant models, the Hall con- 
ductivity must either, if quantized, vanish, or else vary 
continuously and not be quantized. 

The plan of the remainder of this paper is as follows. In 
Section y, we first consider the ground state in a system 
of spinless fermions with p-wave pairing, for the infinite 
plane and for periodic boundary conditions (a torus) . We 
show that a transition occurs between weak and strong- 
pairing phases, which can be distinguished topologically 
in momentum space, or by the number of ground states 
on the torus for even and odd particle number. In Sub- 
sec. [IB, we consider the system in the presence of edges 



and vortices. We argue that there are chiral fermions on 
an edge, and degeneracies due to zero modes on vortices, 
when these occur in the weak-pairing phase. In Subsec. 



H C , we show how the results for ground states can be ex- 
tended to other geometries, such as the sphere. Section O 
as a whole shows that the properties of the weak-pairing 
phase, the ground state degeneracies, chiral edge states 
and degeneracies of vortices agree with those expected in 
the MR phase in the FQHE. The strong-pairing phase has 
the properties expected in the Halperin paired states. In 
Sec. |lll| , we consider the case of spin-triplet pairing, with 
applications to the double-layer FQHE system. There is 
a weak-pairing phase with the properties of the 331 state, 
and also a distinct phase with the properties of the MR 
state. In Sec. IV, we consider spin-singlet d-wave pairing. 
In Subsec. IV A , we argue that the HR state corresponds 
to the transition point between weak- and strong-pairing, 
and so has gapless fermions in the bulk. Then we analyse 
the generic weak-pairing d-wave phase, and argue that it 
corresponds to an abelian FQHE state, with a spin-1/2 
doublet of chiral Dirac fermions on the edge, which has 
also been constructed previously. We also discuss here 
(in Subsec. IV C) general arguments for the existence of 
the edge states and other effects, based on Hall-type con- 
ductivities and induced CS actions in the bulk, for all the 
paired states. In Sec. M, we discuss the effects of disorder 
on all the transitions and phases. The explicit calcula- 
tions of the Hall spin conductivity, in the pure systems, 
are given in the Appendix. 

A brief announcement of our results for pure systems 



was made in Ref. [Q . Some of the results have also been 
found independently by others p2|. 
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II. COMPLEX P-WAVE PAIRING OF SPINLESS 
OR SPIN-POLARIZED FERMIONS 

In this Section, we first set up the BCS effective quasi- 
particle Hamiltonian, and review its solution by Bogoli- 
ubov transformation. We show that this leads to the ex- 
istence of a transition between distinct phases, which we 
label weak- and strong-pairing. They are distinguished 
topologically. The weak-pairing phase is tentatively iden- 
tified with the MR phase because of its behavior in posi- 
tion space. This is extended to the torus (periodic bound- 
ary conditions), where we find three ground states for N 
even, one for N odd, in the weak-pair ing p hase, in agree- 
ment with the MR state. In Subsec. |IIB| , we show that 
the BCS quasiparticles at long wavelengths near the tran- 
sition are relativistic Majorana fermions, and use this to 
analyze the Bogoliubov-de Gennes equations for domain 
walls (edges) and vortices, again arguing that the results 
agre e wi th those obtained for the MR state. In Sub- 
sec. [IC, which may be omitted on a first reading, we 



show how p-wave pairing on a general curved surface can 
be handled mathematically, and that the ground states 
agree with conformal blocks, as expected from MR. 



A. Weak and strong pairing phases 

First we recall the relevant parts of BCS mean field 
theory Q. The effective Hamiltonian for the quasiparti- 
cles is 
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where ^k = £k ~ A* ^^d ^k is the single-particle kinetic 
energy and Ak is the gap function. For the usual fermion 
problems, /x is the chemical potential, but may not have 
this meaning in the FQHE applications. At small k, we 
assume ek — k'^ /2m* where in* is an effective mass, and 
so — /i is simply the small k limit of ^k- For complex 
p-wave pairing, we take Ak to be an eigenfunction of ro- 
tations in k of eigenvalue (two-dimensional angular mo- 
mentum) I — ~\, and thus at small k it generically takes 
the form 



A, 



i\[kx tKy), 



(2) 



where A is a constant. For large k, Ak will go to zero. 
The Ck obey {ck, cj^,} = (5kk'; we work in a square box of 
side L, and consider the role of the boundary conditions 
and more general geometries later. 

The normalized ground state of Kcs has the form 



where |0) is the vacuum containing no fermions. The 
prime on the product indicates that each distinct pair 
k, — k is to be taken once. (We will later consider the 
precise behavior at k = 0.) The functions Uk and Wk sue 
complex and obey |ukp -I- |wkP = 1 to ensure (51|il) = 
1. They are determined by considering the Bogoliubov 
transformation 



«k = "kCk - WkC_k7 

"k = "k4 - ^'kC-k, 



(4) 



so that {ak,ak'} = i^kk' and ak|^) = for all k. By 
insisting that [o^k, ^cff] = -EkCtk for all k, which implies 
that 



Kcs = 22 -^kak"k + const. 



(5) 



with £'k > 0, one obtains (the simplest form of) the 
Bogoliubov-de Gennes (BdG) equations, 



These imply that 
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The functions Wk and Wk are determined only up to an 
overall phase for each k, so they can be multiplied by 
a k-dependent phase, Uk -^ e^'^^uu, Wk -^ e^'^^Vk with- 
out changing any physics. One may adopt a convention 
that one of Wk and Uk is real and positive; in either case 
the other must be odd and of p-wave symmetry under 
rotations. We do not use such a convention explicitly be- 
cause there is no single choice that is convenient for all 
that follows. 

Because of Fermi statistics, which imply (cj^.)^ = 0, we 
can rewrite the ground state (up to a phase factor) as 
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where g^ — fk/uk- Then the wavefunction for the com- 
ponent of the state with A^ fermions {N even) is, up to 
an A^-independent factor. 



*(ri,...,rAr) 



2^/2 (iV/2) 



^E«gn^ 



N/2 



where g{r) is the inverse Fourier transform of gt, 
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and P runs over all permutations of N objects. (For 
fermions with spin, this appears on p. 48 in Ref. [^.) 
The right-hand side of Eq. ( p^ ) is a Pfaffian, which for 
a general N x N matrix with elements My {N even) is 
defined by 



N/2 



FiM 



2W/2(7v/2) 



^sgnP]jMp(2,_i)P(2^), (14) 



or as the square root of the determinant, PfM = 
Vdet M, for M antisymmetric. 

We now consider the form of the solutions to the above 
equations. In the usual BCS problem, the functions At 
and Ek are found self-consistently from the gap equation 
(including Hartree-Fock effects), and /i would be deter- 
mined by specifying the fermion density. However, we are 
not interested in all these details, but in the nature of the 
possible phases and in the transitions between them. We 
expect that the phases can be accessed by changing the 
interactions and other parameters of the problem, but we 
will not address this in detail. In particular, some phases 
may require that the interactions be strong, while BCS 
theory is usually thought of as weak coupling. We will 
nonetheless continue to use the BCS mean field equations 
presented above, as these give the simplest possible view 
of the nature of the phases. 

From Eqs. (p|,|lO|), we see that since -Bk — |Ck| ^ as 
k ^ 0, we will have one of three possibilities for the 
behavior at small k, which will turn out to govern the 
phases. As k — > 0, either (i) ^k > 0, in which case 
|uk| -^ 1, \vk\ — > 0, or (ii) ^k < 0, in which case |uk| — ^ 
0, |fk| — > 1, or (iii) ^k —>■ 0, in which case |uk| and 
|wk| are both nonzero. We will term the first case the 
strong-pairing phase, the second case weak-pairing, while 
the third case, in which the dispersion relation of the 
quasiparticles is gapless, i?k ^ as k ^- 0, is the phase 
transition between the weak and strong pairing phases. 
Thus for fi positive, the system is in the weak pairing 
phase, for fi negative, the strong-pairing phase, and the 
transition is at /i = within our parametrization. The 
reason for these names will be discussed below. 

We now discuss the two phases and the transition in 
more detail. We expect that the large k behavior of ^k 
and Ak that would be produced by solving the full system 
of equations will not be affected by the occurrence of 



the transition which involves the small k's only. Note 
that, at large k, Wk ^ and jUkl -^ 1, which ensures 
in particular that the fermion number, which is governed 
by ^ = I]k"^k, with 



Ck-Ck) = rik 



«k 



(15) 



will converge. Also we assume that i?k does not vanish 
at any other k, which is generically the case in the I — 
— 1 states. Thus within our mean field theory we can 
ignore the dependence of the functions ^k and Ak on the 
distance from the transition, which we can represent by 

In the strong-pairing phase, fi < 0, we have Wk oc 
kx — iky, as k — > 0. Then the leading behavior in 
5k = Wk/uk is (X kx — iky, which is real-analytic in k^ 
and ky. If g^ is real-analytic in a neighborhood of k = 0, 
then ^(r) will fall exponentially for large r, ^(r) ~ e"'"/'"", 
but even if not it will fall rapidly compared with the other 
cases below. Thus we term this phase the strong-pairing 
phase because the pairs in position space are tightly 
bound in this sense. Note that this region is /i < 0, which 
would only be reached for strongly attractive coupling of 
the fermions (we disregard the possibility of other non- 
pairing phases for such couplings). Because nk — l^kP, 
there is little occupation of the small k values in this 
phase. 

In the weak-pairing phase, u^ oc kx + iky for k ^ 0, 
and so ^k c>c l/{kx + iky), which gives 



7(r) oc 1/z 



(16) 



for large r, where z — x + iy. This long tail in ^(r) is 
the reason for the term "weak pairing" . This is exactly 
the behavior of g{r) in the Moore- Read (MR) Pfaffian 
state in the FQHE. In the latter this form, 1/z, is valid 
for all distances. We will therefore try to argue that all 
the universal behavior associated with the MR state is 
generic in the weak-pairing phase, when the theory is 
applied to the paired FQHE states. Notice that in the 
weak-pairing phase, the occupation numbers of the small 
k states approaches 1. Of course, this is also the behavior 
of the Fermi sea. When the attractive coupling is weak, 
one would expect the weak pairing phase. If we imagine 
that only the magnitude of the coupling is varied, then 
when it is small and negative, the BCS weak-coupling 
description is valid, and E^ has a minimum at kp. This 
is not the weak-strong transition. Close to the latter 
transition, _Ek has a minimum at k = 0. As the coupling 
weakens, a point is reached at which the minimum moves 
away from k = 0, and eventually reaches |k| = kp when 
the coupling strength is zero and the transition to the 
Fermi sea (or Fermi liquid phase) takes place. Thus the 
weak-pairing phase does not require that the coupling be 
weak, but is continuously connected to the weak- coupling 
BCS region. 



At the weak-strong transition, /i — 0, we find at small 
k that Ei, = |Ak|, |ukp, |wkP -> 1/2, and gu = {k^ - 
iky)/\\i\. The singular behavior of gk leads to 
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for large r. This behavior is intermediate between those 
of the two phases. Because, like the strong-pairing be- 
havior, it is not a rational function of z, it does not cor- 
respond to a "nice" LLL FQHE wavefunction. Also as 
/i ^- 0, the length scale on which the asympototics of 
the two phases is valid (ro on the strong pairing side) 
diverges. This length should not necessarily be identi- 
fied with the coherence length ^q of the paired state. If 
the latter is viewed as "the size of a pair", it may be 
better associated with the decay of the inverse Fourier 
transform of (ckC_k) = Wk'^k = — Ak/-Ek, which devel- 
ops long-range behavior only at the transition. Also, for 
some purposes the relevant function may be Ak, which 
is always nonsingular at k = 0, and does not vanish at 
the transition. For the purpose of making contact with 
the existing FQHE wavefunctions, ^(r) is the appropriate 
quantity. 

The distinction between the two phases, which we can 
view as requiring the existence of the phase transition, 
does not lie in symmetries, unlike many phase transitions; 
instead, it is topological in nature M. Within the mean 
field treatment, this can be seen in terms of the topology 
of the functions Wk, "^kj or of i5k, Ak. We recall that Wk, 
Wk obey |MkP + |wkP = 1, and that multiplying them both 
by the same phase is irrelevant. These conditions imply 
that they can be viewed as spinor (or "homogeneous") 
coordinates for a 2-sphere S*^, for which two real coordi- 
nates are sufficient for a non-redundant parametrization 
of any neighborhood. As usual, we view u = 1, w = 
as the north pole A/", u — Q, v — 1 as the south pole 
S. We can alternatively parametrize the sphere using 
the unit vector nk = (Re Ak, — Im Ak,Ck)/-Ek, and this 
agrees with the parametrization by Uk, ^k- This is essen- 
tially the "pseudospin" point of view of Anderson [|63[ , 
who described the s-wave BCS state as a spherical do- 
main wall in the pseudospin in k-space, of radius kp. 
Because Uk — > as k ^ oo in any direction, we can add 
a point at infinity in k space and view it also as topo- 
logically a 2-sphere, with the point at infinity at M . The 
functions Mk, t'k thus describe a mapping from S"^ (k- 
space) to S^ (spinor space) , with M always mapping to 
A/". Such maps are classified topologically into equiva- 
lence classes, called homotopy classes, such that maps in 
the same class can be continuously deformed into each 
other. The special cases of maps from 5", n = 1, 2, . . . , 
to any space X define the homotopy groups 7r„(X) (us- 
ing a general method of producing a group structure on 
the equivalence classes with base points, which here are 
M). In our case, 7r2(S'^) = Z, the group of integers. By 
inspection, we find that in the strong-pairing phase, the 



map is topologically trivial. In this phase the map can be 
deformed to the topologically-trivial map that takes all 
k to M ■ In any topologically nontrivial map, as k varies 
over the plane plus point at infinity, u, v (or n) range 
over the whole sphere. Indeed, the number of times that 
a given point ^ Af, such as S, is taken by u, v must be 
at least |to|, for a map in the class labelled by the integer 
m. In our case, in the weak-pairing phase the map passes 
through S at least once, namely when k = 0, and possi- 
bly (most likely) only once, and we can choose to identify 
the class with m = 1. This nontrivial topology in k-space 
associated with non-s-wave weak-pairing ground states in 
two dimensions was pointed out by Volovik g. We also 
note that the p-wave weak-pairing phase map is a topo- 
logically non-trivial texture of the pseudospin that is also 
familiar in physics as a two-dimensional instanton, or a 
skyrmion in a 2 -|- 1 dimensional system, in position space, 
in the 0(3) nonlinear sigma model. Since it is impossible 
to pass smoothly between m = and m = 1, the map 
Mk, Vk must be discontinuous at the transition, which 
is what we found. In fact, |up and |wp tend to 1/2 as 
k — > at the transition, which corresponds to points on 
the equator. The simplest form of such a map is one that 
covers just the northern hemisphere, centered on TV. 

The topological distinction between weak- and strong- 
pairing phases is typical for the fully-gapped complex 
non-s-wave paired states, and not only when the gap 
function is an angular momentum eigenstate. For con- 
trast, note that any s-wave state yields a topologically- 
trivial (to — 0) map. For the s-wave case, it is well known 
that Uk and Vk can be chosen real, and continuously in- 
terpolate between the weak and strong coupling limits 
without a phase transition. Also, in three dimensions, 
the relevant homotopy group for the present spinless case 
is 7r3(S'^) = Z, and so nontrivial states do exist in princi- 
ple, the simplest of which is the Hopf texture in k-space. 
The usual s-wave state is again topologically trivial. Note 
that in the usual BCS model with a separable interaction, 
the gap function vanishes outside a thin shell around the 
Fermi surface, which makes the maps slightly discontin- 
uous. We consider only interactions that are continuous 
in k-space, for which the gap function and the map are 
continuous except at the transition. 

To close this subsection, we consider (as promised ear- 
lier) the effect on the ground states of the boundary con- 
ditions and the total fermion number. That is, we use 
a two-dimensional (Bravais) lattice, and consider a sys- 
tem on a plane with points differing by a lattice vector 
identified, generally described as periodic or generalized 
periodic boundary conditions, or by saying that the sys- 
tem is a torus. 

To obtain a low-energy state in a translationally- 
invariant system when A itself is viewed as a dynam- 
ical parameter, we will assume that A is position- 
independent, and thus we can still use Ak in the quasi- 
particle Hamiltonian in k-space. To be consistent with 



this, the fermions must obey either periodic or antiperidic 
boundary conditions for each of the two primitive direc- 
tions of the lattice (or fundamental cycles on the torus) . 
For a rectangle with x and x + Lx, y and y + Ly iden- 
tified, this means the boundary conditions for the x and 
y directions. These choices of boundary conditions are 
well-known in the description of flux quantization in su- 
perconductors (see e.g. Schrieffer |^). We may imagine 
that either zero or one-half of the flux quantum hc/e 
threads either of the "holes" (fundamental cycles) in the 
torus. The half-flux quanta could be represented either 
by a vector potential, with peridic boundary conditions 
on the fermions in both directions, or by no vector po- 
tential and an antiperiodic boundary condition for each 
direction that wraps around a flux, or by a combina- 
tion of these. The different choices are related by gauge 
transformations. We choose to use boundary conditions 
and no vector potentials, so that A is always position- 
independent. We should be aware that if the gauge field 
(fluxes) are viewed as external, they are fixed as part of 
the definition of the problem, and there will be a sin- 
gle ground state for each of the four possible choices of 
boundary conditions, -|--t-, H — , — h, and , in a nota- 
tion that should be obvious. However, if the gauge field is 
viewed as part of the internal dynamics of the system and 
can fluctuate quantum mechanically (as in highly corre- 
lated systems, including the FQHE, where it is not in- 
terpreted as the ordinary electromagnetic field, and also 
in the usual superconductors where it is) then the four 
sectors we consider correspond to four ground states of 
a single physical system, in a single Hilbert space, albeit 
treated within a mean field approximation. The latter is 
the view we will take. 

For each of the four boundary conditions for the 
fermions, the allowed k values run over the usual sets, 
kx — 2ttVx/ Lx for -t-, iniyx + l/2)/ij, for — , where Vx is 
an integer, and similarly for ky. In particular, k = (0, 0) 
is a member of the set of allowed k only in the case -I-+. 
For a large system, ft and At will be essentially the same 
functions of k for all four boundary conditions, but eval- 
uated only at the allowed values. In the paired ground 
states, k and — k will be either both occupied or both 
unoccupied, to take advantage of the pairing (Ak) term 
in Kqs- When k = is in the set of allowed k, k = 
and — k = cannot both be occupied, because of Fermi 
statistics. However, Ak vanishes at k = 0, so k = will 
be occupied or not depending solely on the sign of ^k=o- 
That is, it will be occupied for /i > (in the weak-pairing 
phase), and unoccupied for /x < (in the strong pairing 
phase), and this is entirely consistent with the limiting 
behavior of ny = |wkp as k ^ in the two phases. At 
the transition, fi — 0, the occupied and unoccupied states 
are degenerate. 

We conclude that in either the weak- or strong-pairing 
phases, there is a total of four ground states, three for 
boundary conditions H — , — h, which are linear com- 



binations of states with even values of the fermion num- 
ber in both phases, but for -|--|- boundary conditions the 
ground state has odd fermion number in the weak-pairing 
phase, even fermion number in the strong-pairing phase, 
because of the occupation of the k = state. In most 
cases, the ground state is as given in Eq. (^), but in the 
weak pairing phase for ++ boundary conditions, it is 
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The ground states specified, whether for N even or odd, 
will have the same energy in the thermodynamic limit 
(not just the same energy density) . Note that if the k = 
state is occupied in the strong-pairing phase, or unoc- 
cupied in the weak-pairing phase, this costs an energy 
Ek=o which we are assuming is nonzero, and all states 
where quasiparticles are created on top of our ground 
states cost a nonzero energy, since we assume that i?k is 
fully gapped in both phases. However, at the transition 
/i = 0, the ground states for ++ with odd and even par- 
ticle number are degenerate, and there is a total of five 
ground states. 

If we now compare with results for the MR state on 
the torus ||32| , [44| , which were derived as exact zero-energy 
ground states of a certain Hamiltonian, then we see that 
the weak-pairing phase for even fermion number agrees 
with the exact result that there are three ground states. 
On the other hand, it was stated in Ref. [Q that there 
are no zero-energy ground states for N odd. Unfortu- 
nately, that result was in error; there is just one such state 
for ++ boundary conditions, which can be constructed 
by analogy with that for the 331 state in Ref. [Q. Before 
turning to the wavefunctions of these states, we also men- 
tion that the behavior found in the present approach in 
the strong-pairing phase agrees with that expected in the 
Halperin point of view ||ll| on the paired states, as Laugh- 
lin states of charge 2 bosons. That point of view predicts 
four ground states for N even, none for N odd. Note that 
in comparing with FQHE states, we factor out the center 
of mass degeneracy which is always the denominator q of 
the filling factor v = p/q (where p and q have no common 
factors) ||64| . The remaining degeneracy in a given phase 
is independent of i^ in the sense that it does not change 
under the operation of vortex attachment, which maps 
a state to another in which 1/v is increased by 1, and 
for generic Hamiltonians this degeneracy is exact only in 
the thermodynamic limit. We note that Greiter et al. 
|p2| claimed that the special Hamiltonian for which the 
MR state is exact should have four ground states on the 
torus for A'' even, though they found only three. They did 
not distinguish the weak- and strong-pairing phases, and 
by assuming that the Halperin point of view is always 
valid, they in effect ascribed the properties of the strong- 
pairing phase to the MR state. In fact, there is a total 
of four ground states in the weak-pairing phase, but the 



fourth is at odd fermion number! They also claimed that 
the statistics would be abelian, even though they consid- 
ered the MR state, and suggested that the prediction of 
nonabelian statistics by Moore and Read ||lj| would hold 
only at some special point. 

If we consider the position-space wavefunctions of the 
paired states on the torus, then for even N in the weak- 
pairing phase we once again find exact agreement of the 
long-distance behavior oi g{ri~rj) with that which holds 
for all distances in the eigenstates of the special Hamil- 
tonian [p2|,E3. Here long-distance means g{r) where r is 
not close to a lattice point (mLxjuLy), m, n integers. 
The long-distance form g r^ 1/z in the plane is replaced 
by an elliptic function (or ratio of Jacobi theta functions) 

11,111 for +-, -+, . The ground state for N odd 

with H — h boundary conditions can be written, similarly 
to one for N even in the 331 case [Q, as 

{N-l)/2 

X n 9++i^Pi2i)-rp(2i+i)) (19) 

Here we can take the torus to have sides L^, L^t in 
the complex plane (Imr > 0, and r — iLy/Lx for the 
rectangle), and 



where A can be complex, and find 
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where di is a Jacobi theta function, and 9[{z\t) = 
d9i{z\T)/dz. g++{r) is periodic because of properties 
of the function 6i mentioned for example in Ref. [[44[ , 
and has a simple pole (7++ ex 1/z as z —^ 0. For the 
ground state with N odd, the non-(complex)-analytic de- 
pendence on y in 9++ cancels. Notice that the unpaired 
fermion with i — P{1) in the terms in ^++ occupies 
the constant, k = single-particle state. When used as 
part of a wavefunction in the LLL, the present function 
is the zero energy state for H — \- boundary conditions on 
the torus for N odd, which was omitted in Ref. iHl. 



B. Majorana fermions, edges and vortices 

In this Subsection, we consider the problems of edges 
and of vortices (which correspond to FQHE quasipar- 
ticles), on which we argue there are chiral fermions 
and zero modes, respectively, in the weak-pairing phase. 
Again, this supports the identification with the MR state. 

We begin by considering in more detail the low-energy 
spectrum near the transition at /i = 0. When fi and k 
are small, we can use 
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which is a relativistic dispersion with |A| playing the role 
of the speed of light. Further, using the same approxi- 
mation, the BdG equations become in position space 



.du ' ^ . 

i-r- = — uu + A I 
at 

.dv 



d_ 
dx 



d 



at 



d_ 
dx 



dy 



M, 



(23) 



1n,y)-, 



(21) 



which is a form of the Dirac equation for a spinor (u, v). 
The BdG equations are compatible with u{r, t) — v{r,t)* , 
and this is related to the existence of only a single fermion 
excitation mode for each k. Thus the quasiparticles are 
their own antiparticles; Dirac fermions with this prop- 
erty are known as Majorana fermions. Near the tran- 
sition, the BCS quasiparticles make up a single Ma- 
jorana fermion quantum field, and at the mean field 
level the critical theory is a single massless Majorana 
fermion. There is a diverging length scale at the tran- 
sition, ^ A/|/i|, and the diverging lengths mentioned in 
the last subsection should all diverge proportionately to 
this, at least within mean field theory. 

Next we wish to consider the behavior near an edge; 
outside the edge the particle number density should go 
to zero. In the Hamiltonian, this can be arranged by 
having a potential V{r) that is large and positive outside 
the edge. In the quasiparticle effective Hamiltonian, this 
can be viewed as making fi large and negative outside 
the edge, and we will use this notation. 

In general, the problem with the edge should be solved 
self-consistently, which involves solving the gap equation 
for A in the presence of the edge. We are interested in the 
generic properties of the solution, and wish to avoid the 
complexities. Accordingly we will consider only a sim- 
plified problem, which is the effective Hamiltonian with 
a given gap function. Since /i becomes negative outside 
the edge, it must change sign near the edge if it is pos- 
itive inside the bulk of the system. But in the bulk at 
least, a change in sign across a line represents a domain 
wall between the weak- and strong-pairing phases, since 
/I = is the point at which the transition occurs in our 
treatment. Thus we are arguing that the weak-pairing 
phase (where fi is positive) must have a domain wall at 
an edge, while the strong-pairing phase need not. We 
will consider a domain wall in the bulk as a model for 
the edge of the weak-pairing phase. In the latter the re- 
gion of strong pairing between the domain wall and the 
actual edge may be extremely narrow and we might say 
there is no well-defined strong-pairing region. But the 
point is the topological distinction between the phases. 
The strong-pairing phase has the same topology as the 
vacuum, and can be continuously connected to it. The 



weak-pairing phase is nontrivial and the generic prop- 
erties of an edge should be captured by a domain waU 

H- 

We consider a straight domain wall parallel to the y- 
axis, with /x(r) = ii{x) small and positive for x > and 
small and negative for x < 0, and also varying slowly 
such that the above long-wavelength approximation can 
be used, with /i now x-dependent. 

We can consider solutions with definite ky, which at 
first we set to zero. Then we have a f + f-dimensional 
Dirac equation. We assume that |/z| — > /uo, a constant, as 
|a;| ^ cx). At £' = there is a normalizable bound state 
solution first obtained by Jackiw and Rebbi [Q. The 
equations are (with A real and > 0) 



Aidv/dx = iJi{x)u, 
Aidu/dx = —ii{x)v. 



(24) 



By putting v = iu, we find the unique normalizable so- 
lution 

1 f^ 
^(x) ex e~*''/''exp-— / fi{x)dx, (25) 

where the phase was inserted to retain v = u* . Solutions 
at finite E should exhibit a gap, as in the bulk in either 
phase. 

At finite ky, the equations become 
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For E = —Aky, these have solutions that are bound to 
the domain wall, and clearly they propagate in one direc- 
tion along the wall. There is only one fermion mode for 
each fcj,, and so we have a chiral Majorana (or Majorana- 
Weyl) fermion field on the domain wall. 

It will be important to consider also a pair of domain 
walls. We consider two walls, lying dX x — Q and x — W, 
with fi>OmO<x<W and /i, < outside. Again 
we assume ky = initially. This time [|66| there are no 
E = Q modes for finite W. Clearly as M^ ^ oo we expect 
to find an E' = mode on either wall, so we expect bound 
solutions for small E when W is large but finite. For non- 
zero E we can replace the pair of first-order equations 
with a single second-order equation for either oi u ± iv, 
(with ky ^ for generality). 
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When /i varies slowly compared with its magnitude /io 
far from the walls, we may study the equations by the 



WKB method. We may view the equation for uztiv as 
a Schrodinger equation with potential 
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If dfji/dx has extrema at x = 0, W , as is reasonable, 
then V-{x) has minima at a; = 0, W ^ but that a,t x — Q 
is deeper than that at x — W. The reverse is true for 
V+(x). As ly ^ oo, there will be a fcj, = solution for 
E'^ which — > exponentially and which corresponds to a 
normalizable eigenfunction for u — iv that is concentrated 
at a; = with negligible weight ai x = W, and similarly 
an eigenfunction for u + iv concentrated at x — W. The 
subtle but important point is that these solutions are 
not independent, because they are related by the original 
first-order system, for any non-zero E. There is only a 
single normalizable solution for the pair u, v for E small 
positive, and another for _E < 0. Consequently, there 
is only a single fermion mode, shared between the two 
domain walls, not one on each. For non-zero fcj,, the 
eigenfunctions for E > become concentrated on one 
wall or the other, depending on the sign of ky. Thus 
the set of low-energy states can be viewed as a single 
non-chiral Majorana fermion theory, with the left-moving 
modes on one wall, the right-moving modes on the other, 
and the ky = mode shared between the two. This 
agrees precisely with the results of Ref. mJ/ on the edge 
states of the MR state on the cylinder. 

We next consider the quasiparticle spectrum in the 
presence of vortices of the order parameter; in two dimen- 
sions vortices are point objects. These necessarily con- 
tain an integer number of half flux quanta in the gauge 
field; without the gauge field the vorticity is quantized 
but the total energy of an isolated vortex diverges loga- 
rithmically. We will not see these effects in the energy 
here because we only consider the quasiparticle excita- 
tion spectrum in the presence of a given gap function and 
gauge field configuration. Up to now we have ignored the 
U(l) gauge field except in discussing the boundary con- 
ditions on the torus and cylinder. It could be the stan- 
dard electromagnetic field in a superconductor, or the 
CS field in the FQHE. We will only consider vortices of 
the minimal flux, namely a half flux quantum, because 
addition of any integer number of flux quanta can be 
viewed, on scales larger than the penetration depth, as 
a gauge transformation, which does not affect the spec- 
trum. Outside the vortex core, which we assume is small, 
the covariant derivative of the gap function must vanish. 
As in the case of the torus, we will choose a gauge in 
which the gap function is single-valued and independent 
of the angle relative to the position of the nearby vortex, 
but the Fermi fields are double-valued on going around 
the vortex. 

The basic idea is to consider a vortex as a small circular 
edge, with vacuum (vanishing density) at the center. Ac- 
cordingly, we expect that nothing interesting happens at 
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sufRciently low energies for vortices in the strong-pairing 
phase. But a vortex in the weak-pairing phase must in- 
clude a concentric circular domain wall to separate the 
vacuum at the center from the weak-pairing phase out- 
side. We now study this using the Majorana fermion 
equations near the transition, assuming that the wall has 
large enough radius; the vortex core (where A vanishes) 
can be taken to have negligible size and the boundary 
condition at r —> is unimportant. With our choice, the 
BdG equations for a single vortex and for E = becomes 
in polar coordinates r, 9 



. jy , i9 id 

dr r 86 

d id 
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u obeys u(r, 6 + 27r) — ~u{r, 9), and similarly for v. We 
can assume fi —^ fj,o > as r —> oo, fj, — > —fio as r — > 0. 
The normalizable solutions have the form 



w=(jz)-i/2/(r), 



(30) 



where fir) is a real function. The equations reduce to 

df/dr = -^i{r)f{r)/A, (31) 

with solution 



/(r) (X exp(- / fi{r')dr'/A) 



(32) 



Thus we find just one normalizable bound state at zero 
energy. Again we expect this to persist as we relax our 
assumptions, as long as the bulk outside the vortex is 
in the weak-pairing phase. We point out that our result 
should be contrasted with the known result for a vortex in 
an s-wave superconductor, which has bound quasiparticlc 
modes at energies that are low in the weak coupling limit, 
but not generally zero as ours are here |67|| . Since we 
mainly work at moderate or strong coupling, analogous 
modes are not important for our purposes. We note that 
a zero mode on a vortex in an A-phase p-wave paired 
state was first found in Ref. |6q| . 

For the case of 2n well-separated vortices, we have not 
obtained analytic solutions for the bound states. How- 
ever, we can give a simple argument. The problem is 
analogous to a double-well potential. We take a set of 
2n E ^ solutions like Eqs. (|30|), (p3) centered at each 
vortex, and use these as a basis set (we must introduce 
additional branch points into each basis state to satisfy 
the boundary conditions at all the other vortices); at 
finite separation there is mixing of the states, and the 
energies split away from zero. Since the solutions to the 
Dirac equation are either zero modes or E, —E pairs, 



we expect to obtain n E > solutions, n E < solu- 
tions. In general, each E > solution of the Dirac or 
BdG equation corresponds to a creation operator, and 
the related E < solution to the adjoint (destruction) 
operator, while an E = solution would correspond to 
a real (or Majorana) fermion operator. In our case, this 
means that there are n modes in which we may create 
fermions, with energies E tending to zero as the separa- 
tion diverges. (A similar picture applies for 2n domain 
walls.) This is in agreement with the results for the spe- 
cial Hamiltonian E9]. This result is crucial for the non- 
abelian statistics we expect in the FQHE case, since by 
occupying the zero modes one obtains a total of 2" de- 
generate states, or 2"~^ for either even or odd fermion 
number TV, when there are 2?! vortices (n > 0); this was 
found for the special Hamiltonian in Ref. |^,Q . 

We may also consider here the edge states of a system 
in the form of a disk of radius i?, by studying a large 
circular domain wall enclosing the weak-pairing phase, 
and strong-pairing phase or vacuum outside. In this 
case, there is no flux enclosed by the wall, and u and 
V are single- valued. One does not find E = states, 
but instead there is a set of chiral fermion modes with 
angular momentum m quantized to half-integral values, 
TO e Z -t- 1/2, and E oc m/R (this fixes the definition of 
m = 0). These are just the modes expected for the chiral 
Majorana fermion on such a domain wall with the ground 
state inside, since an antiperiodic boundary condition is 
natural for the ground state sector. If a half-flux quan- 
tum is added at the center of the disk, the quantization 
is TO, S Z, and this extends the result for the zero mode 
771 = of a single vortex. These results agree with the 
results of Ref ^^JU^ for a disk of the MR state. 

We also note that the form of the modes near a vor- 
tex, containing z~^" or its conjugate, is similar to the 
form of the fermion zero mode functions found in Ref. 
[Q, once the factors associated with the charge sector 
are removed, though the factors /(r) are not. However, 
the most appropriate comparison to make is that be- 
tween the many-fermion wavefunctions, as we already 
made for the ground states on the plane and torus. We 
will not consider this further here for the vortex or the 
edge (or domain wall) states, though we expect that these 
should correspond at long distances to those found in Ref. 
P3| , p4|J42[ , as for the ground states. However, we are able 
to find the ground states on other geometries, namely the 
sphere and Riemann surfaces of genus (number of han- 
dles) greater than one (the sphere is genus zero, the torus 
genus one). We consider this briefly in Subsection II C 
below. 



C. Other geometries and conformal field theory 

In this Subsection, we briefly introduce some general 
connections of pairing theory for p-wave states to rela- 
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tivistic fermions, which enables us to discuss geometries 
other than the plane and torus, such as the sphere, and 
to make more explicit connections with conformal field 
theory ideas iQ. This Subsection can be omitted on 
a first reading, but some of the formalism is mentioned 
again later. 

In the preceding Subsection, we used the fact that the 
BdG equations at long wavelengths become the Dirac 
equation, with a reality condition so that the Fermi field 
is a Majorana fermion. We also mentioned the coupling 
of the gap function and Fermi fields to a U(l) vector 
potential (which in the FQHE context would be the CS 
vector potential), which is of a standard form. But the 
interpretation of the fermion as Majorana would seem 
to raise a problem, because for a single Majorana there 
is no continuous symmetry of the Yang-Mills type, and 
so apparently no way to minimally couple it to a vector 
potential. We will see that there is nonetheless a natural 
way to incorporate the vector potential and still give an 
interpretation in terms of the Dirac equation, and this 
will also enable us to discuss the ground states on curved 
surfaces. 

The most general form for the p-wave gap function 
in Fourier space, retaining once again only the long- 
wavelength part, can be written 



ZAIc — IXx.f^x 



L '\t 1^11 < 



(33) 



Here ^x,y are two complex coefficients, or equivalently 
four real numbers, which we will arrange into a 2 x 2 ma- 
trix e. In position space, the k can be replaced by — iV. 
Then in a general coordinate system x', with correspond- 
ing partial derivatives i9i, ii = x, y), the BdG equations 
become 

{dt + izw^Efe,)^ + e^'^aaid, + \^iJ^^ll^M -K ^/3mV^ = 0, 

(34) 

where we use a spinor i/) — {v,u), and m — fi in previous 
notation. The indices a, b, c take the values x, y, and 
the matrices are a^, = a^, cty = (Jy, P = a^', we use the 
summation convention. Here we have also reinstated the 
vector potential A^ — w^^ /2 (where fi — t, x, y), using 
the matrix 



E. 



(35) 



The equation is therefore invariant under U(l) gauge 



transformations ?/; 



oiA'^S. 



"i/'j and a correspond- 



ing transformation of w^, with a real scalar parameter 

If we multiply through by /3, then we obtain the more 
covariant form of the Dirac equation, 

e'^"7a(9^ + ^<Efc,)^ + tm^j = 0, (36) 



with e'^* = Sfj^t in our case. This has the form of the 
general Dirac equation suitable for use in general coor- 
dinate systems on general curved spaces or spacetimes 
[ |69[ . A similar form was obtained in Ref. 0]. In gen- 
eral, -0 is a Dirac spinor (with two components in the 
2 and 2+1 cases of interest here); /i = 1, 2, ... d is a 
spacetime index, while a, 6, c = 1, 2, . . . , d is an inter- 
nal "local Lorentz" index; the vielbein e is a tensor with 
indices as shown; 7^ are a set of Dirac matrices satis- 
fying {7a, 76} — 'irjab, where 77 is the Minkowski or Eu- 
clidean metric, and Hab = ^i[7a,7b] are the generators of 
SO(d— 1, 1) Lorentz transformations (or simply SO((i) ro- 
tations, in the Euclidean case); the spin connection w is a 
tensor field with one spacetime and two internal Lorentz 
indices, and is antisymmetric in the latter. Spacetime 
indices are raised and lowered using g'^'^ and 5^1/, while 
internal Lorentz indices are raised and lowered using 77"'' 
and rjab- This form of the Dirac equation is covariant un- 
der coordinate transformations (diffcomorphisms) , under 
which the spinor is viewed as transforming as a scalar 
function of position, and e and uj as tensors. It is also co- 
variant under SO((i— 1, 1) (or SO{d)) local Lorentz trans- 
formations, which act like gauge transformations, with 
■0 transforming in the spinor representation, the vielbein 
transforming as a vector in the a index, and uo transform- 
ing inhomogeneously as a nonabelian vector potential or 
connection for the gauge transformations. This formal- 
ism (also known, in four dimensions, as the vierbein or 
tetrad formalism) can be used to reformulate, for exam- 
ple, general relativity in a form equivalent to the usual 
one involving Christoffel symbols; this involves imposing 
relations 






(37) 



On the other hand, it is the only known way to cou- 
ple Dirac fields to curved spacetime; more details can be 
found in Ref. |69|. In our case, we have a distinguished 
choice of time coordinate, we consider only the restricted 
form with e^* — 5^n and require covariance under only 
the S0(2)?^U(1) subgroup that describes "internal spa- 
tial rotations" . Then the Dirac equation becomes pre- 
cisely the BdG equation. This relation with the vielbein 
formalism suggests that the vector potential will play a 
natural role when we consider pairing of nonrelativistic 
fermions on a curved surface. 

The problem of spinless p-wave pairing of nonrelativis- 
tic fermions moving on a general curved manifold should 
be formulated as follows. The manifold has a metric and 
a corresponding Riemann curvature tensor, which for two 
dimensions reduces to a curvature scalar. We will con- 
sider only the case in which this curvature is constant 
on the manifold, and we will also introduce a U(l) or 
SO (2) gauge potential. For a large system, the radius 
of curvature of the manifold is large, and locally the so- 
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lution to the gap equation should resemble that in flat 
space, which we assume is of the I = —1 form. In order 
to minimise the energy density, we expect that the gap 
function should be as constant as possible. In more in- 
variant language, this means that the vielbein should be 
covariantly constant (note that the covariant derivative 
of the vielbein must be covariant under both coordinate 
and U(l) gauge transformations), and there should be 
no vortices. Since the Riemannian geometry (the metric 
and the Levi-Civita or metric connection) of the manifold 
are assumed given, this condition relates the spin connec- 
tion Lo to derivatives of the vielbein, a result analogous to 
the usual requirement (for s-wave pairing, in flat space) 
that the vector potential be the gradient of the phase 
of the gap function. This covariant-constancy condition 
also appears when formulating general relativity p9\ . 

This condition can be satisfied globally (we assume the 
surface on which we are working is compact) only if the 
field strength in the SO (2) vector potential is related to 
the Riemann curvature of the manifold. The integral of 
the latter over the surface, divided by 47r, is a topologi- 
cal invariant, the Euler invariant, equal to 2(1 — 17) for ^ 
Riemann surface of genus g (one with g handles). In our 
usual units for flux quanta, the number of flux quanta in 
the S0(2) or U(l) vector potential must he g — 1. Oth- 
erwise, we will have vortices somewhere on the surface, 
at which A goes to zero. In particular, for the sphere, 
this agrees with the familiar fact for the MR state that 
the number of flux iV^ seen by the underlying particles 
is one less than in the Laughlin state at the same fill- 
ing factor, so the composite fermions see a net flux of 
— 1. Physically, the nonzero angular momentum of the 
pairs causes them to see the curvature of the manifold 
on which they move as a gauge fleld, the fleld strength 
of which is cancelled (locally, not just globally) by the 
imposed U(l) (i.e. S0(2)) gauge fleld, so that a uniform 
condensate is possible, much like the condition of van- 
ishing fleld strength for uniform s-wave condensates in 
ordinary superconductors. 

It can be shown that the long- wavelength wavefunction 
involves the inverse of (or Green's function for) part of 
the (covariant) Dirac operator we have discussed, namely 
the part A^, where A is the part of the Dirac operator, 
including the vector potential, that acts on ct and maps it 
to c (like the earlier gap function) in the Dirac equation. 
A^ contains derivatives like d/dz in local coordinates. 
On the sphere, in stereographic coordinates the Green's 
function is known to be essentially l/{zi — Zj) for parti- 
cles i and j. For any surface, this description in terms of 
inverting the massless Dirac operator is identical to the 
problem of flnding the correlators of two-dimensional chi- 
ral Majorana fermions (in Euclidean spacetime), and so it 
is not surprising that this agrees with the conformal block 
for N two-dimensional massless fermions on the sphere in 
conformal fleld theory. We note that the paired ground 
state on the sphere can be described in angular momen- 



tum space, in terms of single-particle angular-momentum 
eigenstates with eigenvalues j, m (and j = 1/2, 3/2, . . . , 
due to the single flux quantum), as BCS pairing of j, m 
with j, —m; antisymmetry and vanishing total angular 
momentum for each pair require that the j's be half-odd- 
integral, as they are for p-wave. 

We may also consider the cases of Riemann surfaces 
of genus greater than one. Here the explicit functions, 
which again are built out of 1/At, are more difflcult to 
flnd, but certainly exist and describe the MR state on 
these surfaces. The required number of flux seen by the 
fermions is g — 1. The only aspect we wish to discuss 
further here is the number of distinct ground states for 
g > 1. When handles are present, the vector potential 
is determined only up to addition of a pure gauge piece 
describing holonomy around the 2g fundamental cycles 
of the surface. The holonomy, or phase picked up when a 
fermion is parallel-transported around a cycle, can only 
be ±1, since it comes from the double cover of SO (2) 
by Spin(2), the group which possesses the spinor rep- 
resentation. This effect, which is a restatement of flux 
quantization, agrees with and generalizes the discussion 
of boundary condition sectors for the torus. There is thus 
a set of 2^^ possible boundary condition sectors, which 
in the present differential geometry set-up are known as 
spin structures. The spin structures on a genus g surface 
can be divided into two sets, known as the even and odd 
spin structures. The difference between these, for our 
purposes, is that the odd spin structures possess a sin- 
gle zero mode for the Dirac operator, and the even spin 
structures possess none. Then the BCS ground states in 
the weak-pairing phase will include one fermion occupy- 
ing the zero mode when one exists, and since the other 
fermions are all paired, we conclude that the odd spin 
structures give rise to ground states with N odd, and 
the even spin structures to N even, and these ground 
states will be degenerate in the thermodynamic limit. It 
is known that there are 2^"-'^ (2^ -I- 1) even spin structures, 
and 2^^^ (2^ — 1) odd spin structures, so these formulas 
give the number of ground states for TV even and odd, 
for all g > 0. These numbers (and the long-distance 
wavefunctions) agree with the conformal blocks for a cor- 
relator on the genus g surface with N Majorana flelds 
inserted. All of this is in beautiful agreement with the 
CFT picture of Ref. ||l|l. We note that the U(l) charge 
sector which is present in the FQHE states gives another 
factor q^ in the degeneracy for flUing factor 1/ = p/q, 
in the thermodynamic limit [|70|, which generalizes the 
center-of-mass degeneracy q of the torus pM . 



III. SPIN-TRIPLET COMPLEX P-WAVE 
PAIRING 

In this Section we consider spin-triplet p-wave pairing. 
Since the general classification of such states is compli- 
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cated (compare the three-dimensional version in Ref. Q), 
we concentrate on a particular case directly related to 
the FQHE. The FQHE system we have in mind is the 
double- layer system at v = 1/2 [Q. This is assumed 
to be spin-polarized, but the layer index of the electrons 
plays the role of a spin, which we refer to as isospin, 
to avoid confusion with the pseudospin discussed in the 
previous Section. The iSz-values of the isospin will be 
denoted t, i for the two layers. 

In the double-layer FQHE system, as in the other sys- 
tems we discuss, we go to a CS fermion representation 
by using layer-independent fluxes attached to the CS 
fermions. Because interactions between electrons in the 
same and in different layers are different (though the two 
layers are on an equal footing), the Hamiltonian will not 
have SU(2) symmetry. However, in the absence of a tun- 
neling term and of interactions that transfer electrons be- 
tween layers, the number N^ — N^ is conserved, and since 
this quantity is twice the total Sz of the isospin, there is 
a U(l) = S0(2) symmetry that rotates the isospin about 
the z axis. Also interchange of the two layers (or re- 
flection in the plane midway between the two layers) is 
a Z2 symmetry. Together these make up an 0(2) sym- 
metry. We also consider the effect of a tunneling term 
—tax for each particle; t is the tunneling amplitude and 
ax, etc, denote the Pauli matrices. Nonzero t breaks the 
symmetry to the Z2 of layer exchange. 

The FQHE system &X, v — 1/2 has a possible ground 
state which is a 331 state, which can be viewed as com- 
plex p-wave pairing of composite fermions EJ . The pair- 
ing and the effective quasiparticle Hamiltonian are best 
considered in terms of isospin states which are eigenstates 
of ax, namely e = (t + i)/\/2, o = (t - i)/v^, which 
are respectively even, odd under the Z2. As in the ear- 
lier work on this problem, we assume that the Z2 sym- 
metry, and for t = the 0(2) symmetry, are not broken 
spontaneously. Then symmetry dictates that the effec- 
tive quasiparticle Hamiltonian has the form 



k L 

+^kecLcLke) + (^k + i)4oCko 

+ 2 (AkoC-koCko + Ako4„C^^ko 



(38) 



We have taken the same kinetic term ^k — fc^/2TO* — /i 
for both e and o since a difference here is unimportant 
(and forbidden by symmetry when t = 0). For t = 0, 
Ake = Ako, and in general we assume both have p-wave 
symmetry, with 



Ake — ^e\f^x ItZy 



(39) 



ting 2i between e and o (such as an exchange enhance- 
ment). 

We see that the unbroken Z2 symmetry has led to de- 
coupled e and o Hamiltonians. These are the same as for 
the spinless p-wave case. Consequently, we see that sepa- 
rate transitions from weak to strong pairing are possible 
when i ^ 0. For i = 0, these coincide. The pairing func- 
tion g{y) in the wavefunction is now a four-component 
object because of the isospin variables. We write it as 
a vector in the tensor product space of the two spinors. 
When at least one of the two components is in its weak- 
pairing phase, the pairing function at long distances has 
the form 



cos(^ — 7r/4)eiej -I- sin(0 — 7r/4)oiOj 



(40) 



This is the form that was assumed for all distances in 
Rcfs. ||J7^,||. For i = 0, we would put 6* = 0, in which 
case it reduces to 



1 j-ij 



u^i 



at small k, and similarly for Ako- We have also neglected 
the possibility of many-body renormalization of the split- 



(41) 



which is the form in the 331 state. As t increases, neglect- 
ing the likely change in Ake, Ako momentarily, a point 
is reached at which the o spins have an effective chemical 
potential 11 — t — 0, and undergo a transition to strong 
pairing. Then the long distance behavior is eiCj /{zi — Zj), 
so 9 — 7r/4 at the transition and remains at that value 
thereafter. The resulting phase is expected to have the 
statistics properties of the MR state, unaffected by the 
strong-pairing 00 pairs present in the ground state. 

When /i is decreased, there will be a transition from 
the MR phase to strong pairing in both components when 
fi + t = 0. Thus we obtain the phase diagram shown in 
Fig. H We have also included labels of the analogous 
phases in He'^ [B[. In He^, the roles of e, o are played 
by Tj ii and that of t is played by the Zeeman splitting 
due to a field h along the z direction. In He"^, there 
is full SU(2) symmetry of spin rotations when h = 0, 
that is broken spontaneously in any spin-triplet phase, 
but this distinction is unimportant here. The state for 
t = has the structure of the ABM state or A phase, 
adapted to two dimensions, while for t ^ 0, we expect 
that self-consistent solution of the gap equation would 
give Ae > Ao, and this state has the structure of the A2 
phase of He^. As t increases, a point may be reached at 
which Ako = for all k, which gives the Al phase (we 
ignore intermediate possibilities in which Ako vanishes 
only in some region of k space). For /i — f > 0, which 
would be the case in He^, the Al phase is a distinct phase, 
which would have a Fermi surface for o spins. On the 
other hand, for /x — i < 0, no excitations become gapless 
at the point where Ako vanishes, and the change is merely 
the disappearance of 00 pairs from the ground state, so 
this is not a true phase transition; this is indicated by 
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FIG. 1. Schematic phase diagram for the p-wave phases, 
as discussed in the text. The A-phase with unbroken U(l) 
symmetry appears as the vertical axis t = 0, with the region 
/i > being the 331 phase. Similarly, the Fermi-liquid phase 
in which pairing disappears is identified with the line ^ = Ef, 
since that is the value of /j, there at fixed density, neglecting 
Hartree-Fock corrections. 

the dashed line in Fig. |[ The position of this boundary, 
if it occurs at all, is very uncertain. We simply wish 
to emphasize that the boundaries between A, A2 and 
Al in general do not coincide with the weak- to strong- 
pairing transitions. However, for the wavefunctions that 
are of the form Pig with g given by Eq. (Eoh for all 
distances, the transition at 6* = n/A to the MR state can 
also be considered as the A2-A1 transition. Wo discuss 
this further below. 

For t — 0, the quasiparticle excitations near k = and 
/.t ~ form a Dirac relativistic fermion spectrum, which 
has two distinct but degenerate particle and antiparticlc 
excitations. These are eigenstates of the U(l) symme- 
try of eigenvalues -fl, —1, respectively, so the particle 
is associated with one layer, the antiparticle with the 
other. A Dirac fermion field is equivalent to two Majo- 
rana fields, which can be thought of roughly as its real 
and imaginary components. For t ^ 0, these Majorana 
fermions have different masses (and different velocities if 
Ae ^ Ao), and the transitions occur when one or other 
mass changes sign. 

The ground-state degeneracies, edge state and vortex 
(quasiparticle) properties in the different phases can now 
be read off from the above and the results for the spin- 
less p-wave case. For /i > and t — 0, we find the 
same results as for the 331 state [|2|j4|]. The FQHE 
state is abelian and can be described in Coulomb plasma 
language because of the Cauchy determinant identity, 
as explained in Ref. [E2lE4|. The equivalence of the re- 
sults for the neutral edge excitations, or for the quasihole 
states, in terms of bosonic fields, or the two-component 
plasma mapping, and Dirac fermion fields, corresponds 
to bosonization p2p3] and will not be repeated here in 



full. We will mention only an instructive example, con- 
sisting of a pair of vortices in the weak pairing phase 
/i > 0. Given that there is a single E — mode for 
the pair, which can here be occupied by either type of 
fermion, then there are four states. Since the fermions in 
the "f-J, basis have Sz quantum numbers ±1/2, the four 
states have Sz values 1/2, 0, 0, —1/2 (half of these states 
have odd total fermion number N). These states must 
be interpreted as saying that each of these elementary 
vortices carries Sz of ±1/4, so there is a fractionaliza- 
tion of the Sz quantum number. For N even we have 
only two states, both with Sz — 0, if there are no other 
excitations of the system, due to global selection rules 
related to the total quantum numbers. This agrees with 
the two-component plasma description for the 331 states, 
on including the charge degree of freedom in the incom- 
pressible FHQE system [ p2|j4^ ] . In the latter formulation, 
the fractional Sz is analogous to the fractional charge of 
the Laughlin quasiparticles. 

For t ^ 0, the U(l) symmetry is lost, and the quantum 
field theories, whether the massive theory in the bulk 
or the chiral theory on the edge or at the vortices, in 
the phase labelled "weak-pairing abelian" must be de- 
scribed as two Majorana fermions. However, the count- 
ing of the edge excitations, or of the vortex states just 
discussed, will be the same (though unimportant degen- 
eracies among edge excitations, that previously were due 
to the U(l) symmetry, may be lost), and the universal 
statistics properties in this phase, which are abelian, are 
the same as in the U(l) symmetric or Dirac fermion case. 
As already mentioned, when the transition to strong- 
pairing occurs for the o spins, the system enters the MR 
nonabelian phase, and when the e spins are also in the 
strong-pairing phase, the FQHE system has the rather 
trivial abelian statistics of the Laughlin state of charge 2 
bosons. 

Ideas of Ho [Q involving a two-body pseudopotential 
Hamiltonian, the ground state of which interpolates be- 
tween the 331 state at 6* = and the MR state at = 7r/4 
were discussed critically in Ref. |Q. The model Hamil- 
tonian contains a parameter corresponding to 6 which 
directly determines the 6 that describes the ground state 
wavfunction. Except at 9 = 7r/4, the results of the model 
agree with the above discussion of the abelian phase, 
though we do not now believe that the U(l) symmetry 
is "reappearing in the low-energy properties" p4[| . At 
6 = 7r/4, Ho's model is pathological, and the upshot of 
the discussion was that if the three-body interaction of 
Greiter et al. is also added as a isospin-independent in- 
teraction, then this pathology of Ho's model is removed, 
and the ground state is nondegenerate (on the sphere, 
without quasiholcs) for all 6 between and tt/4. It was 
noted that, even for this model, there are still peculiari- 
ties of the 6 = 7r/4 point. In fact, there are unexpected 
degeneracies, larger than those of the 331 or < 9 < 7r/4 
states, in the edge, quasihole and toroidal ground states 
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at 6 = 7r/4. These arise because the zero modes can be 
occupied by either e or o fermions, even though there are 
no OS in the ground state. In particular, on the torus 
in the ++ sector, there are degenerate states that dif- 
fer only in the presence or absence of a k = o fermion 
(these extra degeneracies in this model at this point were 
overlooked in Rcf. Um , but can be obtained by the same 
methods used there). This clearly suggests that the dis- 
persion relation Ei^^o for the o Majorana fermions is gap- 
less at this point, and seems to confirm that this model 
at 9 = IT /A is actually at the transition point between the 
weak-pairing abelian and MR nonahelian phases. This 
is consistent with the result of the analysis here that 
when the leading long-distance part of g is described by 
9 = 7r/4, the system is either at the transition or in the 
MR phase. This is certainly not the case for 9 < 7r/4, as 
indicated by the degeneracies found for the Ho Hamilto- 
nian plus three-body in Ref. Q]. However, the fact that 
the ground state of the model at 6* = 7r/4 contains no o 
fermion pairs at all, suggests that this point is at the A2- 
Al boundary as well as at the weak-MR transition. In 
the quasiparticle effective Hamiltonian, this would cor- 
respond to vanishing Aq as well as fi, and therefore we 
would expect the dispersion relation for the o fermions to 
be -Ek oc |kp. Clearly this is nongeneric behavior. When 
tunneling t is also included, the value of 9 in the ground 
state cannot be read off the Hamiltonian in general, but 
we expect that the weak to MR transition, which should 
now be generic, is pushed to another value of the param- 
eter in the Hamiltonian, so the region in the MR phase 
with 9 — 7r/4 widens, while the A2-A1 boundary is still at 
the value corresponding to 9 = 7r/4 originally, where the 
ground state is known exactly p4[ , and again contains 
no o fermions. 



IV. COMPLEX D-WAVE SPIN-SINGLET 
PAIRING 

In this Section we consider I — —2 complex d-wave 
pairing of fermions, and more generally dx2_y2 + idxy, 
which are necessarily spin-singlet. This has been con- 
sidered recently [p9|j5q - pq| . We argue that the Haldane- 
Rezayi (HR) state [p_3[\ which has this symmetry, is at the 
transition from weak to strong pairing. The weak pairing 
phase, like the strong, is abelian; we work out its prop- 
erties and its lattice description from the pairing point 
of view, and identify the universality class as one which 
has been obtained before by various methods, including 
a trial wavefunction approach by Jain [^2|J7|-Q. We 
relate these properties to the Hall conductivity for spin, 
which we calculate explicitly in an Appendix. 



A. Weak- and strong-pairing phases, and the 
Haldane-Rezayi FQHE state 

The basic structure of the problem is once again similar 
to the spinless p-wave case, or (since it is spin-singlet) 
to the original BCS treatment, except for being d-wave. 
The quasiparticle effective Hamiltonian is 



K, 



(CkcLck<T 



A^c_kiCkT + Ak4-fcLk J , (42) 



where Ak — A{kx — ikyY at small k. The structure of 
the solution is similar to the spinless case. However, the 
dispersion relation is now 



E^,= 



/(fcV2^^^^~77"7S^, 



(43) 



so at the transition point [ji = 0), E\^ = fc^(|Ap 
(2to*)^^)^/^. For ^ nonzero, 



E^ - \^i\ 



2m 



rsgn/ii. 



(44) 



at small k, which implies there is a minimum at nonzero 
k in the weak-pairing phase, /i > 0. 

The position-space wavefunction for N particles, of 
which N/2 have spin "f, A^/2 have j, has the form 



*cxdet5(ri| -r^x), 



(45) 



where g is the inverse Fourier transform of gk = Wk/uk 
(this is equivalent to a result in Ref. [g|, p. 48). In the 
strong-pairing phase, fk/wk ~ {k^ — iky)^, and ^(r) falls 
rapidly with r. In the weak-pairing phase, f k/uk ^ (kx + 
iky)~'^ , we find 



g(r) ex z/z 



(46) 



for large r. Thus \g\ ^ constant, and g is very long-range. 
At the critical point, t^k/wk ^ |kp/(fca; + ikyY (with a 
coefficient that depends on A and m.*, unlike the p-wave 
case) and 



g{v) ex 1/z^ 



(47) 



This is the same behavior as in the Haldane-Rezayi (HR) 
state |13| , when the latter is interpreted in terms of pair- 
ing of composite fermions [H . Therefore we suggest that 
the HR state is precisely at the weak-strong pairing tran- 
sition point, and has gapless excitations in the hulk. 

Further evidence for the criticality of the HR state 
comes from the ground states on the torus. For the quasi- 
particle effective Hamiltonian KcSi the presence of two 
spin states means that the k = states, which occur only 
for -|--|- boundary conditions (see Section O), can be unoc- 
cupied in the strong-pairing phase, and doubly-occupied 
in the weak-pairing phase. Thus there is a total of four 
ground states, all with N even, and none with N odd, in 
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both phases. However, at the critical point, -Ek=o — 0, 
and the k = state can be occupied zero, one, or two 
times, with no energy penalty. Hence for ++ boundary 
conditions, there are two ground states for ISi even, and 
two for N odd. The latter pair, in which k = is singly 
occupied, form a spin- 1/2 doublet. Thus (on including 
the three ground states in the other sectors) there is a to- 
tal of five ground states for N even, two for N odd. This 
is exactly what was found for the ground states of the 
hollow-core model of HR |13| , (for which the HR state is 
the exact ground state on the sphere), both numerically 
[ f75| , p[ and analytically Q. The long distance behav- 
ior of the wavefunctions implied by the present approach 
agrees with that in Ref. ^J], as in earlier cases. Also, 
we cite the energy spectrum of the hollow-core model, 
which was obtained numerically for N — 8 particles on 
the sphere in Ref. [Q. No clear gap can be seen in the 
spectrum. In view of these results, other numerical work 
on this model should also be reconsidered. Analytical 
results on zero-energy "edge" and "quasihole" states of 
the hollow core model [^ 44|, remain valid, but the ear- 
lier interpretation assumed a fully-gapped bulk ground 
state, and so if the bulk is gapless, the questions about 
the conformal field theory pictures of the bulk wavefunc- 
tions and the fermionic edge excitations |14| , pd| , p^p5| ] are 
presumably moot. It is actually quite interesting that the 
hollow-core model (like the Ho plus three-body model in 
Section [II) is critical. In this it resembles also certain 
other special Hamiltonians for which the exact ground 
states are known pGlQ; the latter cases involve pairing 
of composite bosons. 



B. Structure of the weak-pairing abelian phase 

It is now of interest to find the properties of the d- 
wave weak-pairing phase (the strong pairing phase has 
the same rather trivial properties as the others discussed 
previously). This phase has been discussed recently [p7[ . 
Here we wish to consider its application as a paired state 
in the FQHE. There are differences in the symmetry 
here compared with Ref. |£j, which necessitate a cer- 
tain amount of discussion. As we will see, the method 
of analysis of the edge and vortex states used previously 
does not seem satisfactory in the present case. There- 
fore, it seems necessary to use a more devious approach, 
which we now describe. 

It is not difficult to see that, in the weak-pairing phase, 
the map from k to Wk, Vk is topologically non-trivial and 
has m — 2, that is, it wraps around the sphere twice. 
Because of this and its d-wave {I = —2) symmetry, it has 
double zeroes at k = and k = oo. While such behavior 
at oo can be regarded as fixed by requirements of conver- 
gence and finite particle number, as mentioned earlier, 
that at k = is non-generic from the topological point 
of view. More generically, the map could pass over the 



south pole S in Anderson pseudospin space at two differ- 
ent k values, but this would require that the rotational 
symmetry be broken. 

It will be useful to analyze such generic behavior in 
order to find the properties even of the d-wave case, as we 
will argue below. A convenient way to break rotational 
symmetry is to introduce an s-wave component Aks of 
the gap function Ak, in addition to the d-wave part Akd 
(s-wave is the simplest choice, and the gap function must 
remain even in k, to retain the spin-singlet ground state). 
Then at small k the behavior is 



Ak ~ A, -I- A{k^ - ikyf 



(48) 



and As and A are both complex. In this case, when Aks 
is a sufficiently small perturbation on the weak-pairing 
d-wave phase, the effect is that _Ek has a minimum at 
just two nonzero values of k, of the form ±ko by sym- 
metry. This does not change the topology of the Wk, 
fk, and since the system is fully gapped as far as the 
quasiparticles are concerned, it should not change any 
physical properties. However, if Aks is large and domi- 
nates Akd for all k, then the system is essentially in an 
s-wave state which must be topologically trivial. This is 
the same phase as the strong-pairing d-wave phase, even 
though the rotational symmetry is different; these limits 
can be connected as Aks, Akd, and /z are varied, with- 
out crossing another transition. Therefore there must be 
a phase transition as |Aks/Akrf| varies at fixed /x > 0, 
similar to those discussed above. In this case, the disper- 
sion relation E\^ vanishes linearly at two points ±ko at 
the transition, and the map from k to Uk, I'k is discon- 
tinuous at these two points. These will be points where 
Ck = k'^ /2m* — fjL — Q, as well as Ak — 0, and so occur at 
some // > 0. As /i — + 0, ko — + 0, and these points coalesce 
to give the previous discussion in which Aks = for all 
k. 

The two points ±ko at which _Bk has a conical form 
give a spectrum similar to a spin-1/2 doublet of Dirac 
fermions. This is similar to behavior well-known in 
other condensed-matter models, including fermions on 
a lattice in a magnetic field, and d-wave pairing with a 
dj.2_y2 structure (perhaps induced by a square lattice) 
pTj . By concentrating on the degrees of freedom near 
these points, and shifting them in k-space to the origin 
(which produces oscillating factors in real space correla- 
tion functions of the fermions), the fermion excitations 
can be mapped onto a complex (i.e. Dirac, not Majo- 
rana) spin-1/2 doublet of relativistic fermions for each 
pair of such points in k space. Near the transition point, 
or by varying parameters in the other models mentioned 
(the dxy part of the gap function in the second model), 
the Dirac fermions gain a mass. If we now apply to this 
a similar analysis for edges to that we used previously 
when the minima were at k = 0, then we find that the 
edge excitations form a spin-1/2 doublet of chiral Dirac 
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(or Weyl) fermion fields |£^, and we can give a similar 
analysis for the vortices. 

For our purposes, we are interested in unbroken rota- 
tional symmetry in the bulk, so the preceding remarks do 
not seem to apply directly. However, if we analyze the 
edge excitations using the method of previous sections 
applied to the d-wave model, by simply considering the 
effect of ^ changing sign on a domain wall, we do not 
find any edge modes. But we are suspicious of this result 
because of the nongeneric form of Uk, Wk, and previously 
we were relying on the assumption that the results are ro- 
bust because of the topological nature of the phenomena. 
In the present case, the results should be the same as if 
we examine a domain wall caused by varying Aks so as 
to cross the transition, because the phases are the same, 
and therefore the edge state and vortex properties should 
be those of the Dirac fermion doublet. We expect that 
what happens, even if the bulk phases are rotationally in- 
variant, is that the edge or domain wall breaks rotational 
symmetry, and induces a splitting of the zeroes of Ak, in 
so far as this function in k space is meaningful. With the 
symmetry broken, the previous analysis can be applied. 
Of course, what should be done in all cases is a full so- 
lution of the BdG equations with the reduced symmetry, 
and of the gap equation self-consistently. This is clearly 
difhcult, though it has sometimes been attempted, and a 
deeper analysis that explains why arguments of the type 
we have given yield the correct results would be prefer- 
able. We will attempt to give such an argument below. 
We note that, if we consider perturbations analogous to 
the s-wave component in the cases in earlier sections, we 
do not find any change in those results. 

A feature of the Dirac-like nature of the fermion edge 
spectrum is that there are two doublets of fermion modes, 
particles and antiparticles. These correspond to the two 
points ±ko in k space from which they arose. This can be 
described by saying that there is a U(l) quantum number 
or "charge" , which we will call Af , and that the particle 
and antiparticle carry opposite values of M. It may be 
that the additional U(l) symmetry, which must not be 
confused either with any part of the SU(2) of spin, or 
with the U(l) of underlying particle- number conserva- 
tion, which is already broken spontaneously by the pair- 
ing, is in fact broken by the dynamics, since there seems 
to be no symmetry to protect it. However, as in the case 
of the 331 state plus tunneling, even when the symme- 
try is broken, the counting of edge excitations, and the 
statistics properties, etc, should be unchanged. 

If we continue to assume the U(l) symmetry exists, for 
the sake of the latter analysis, then the doublet of chiral 
Dirac fermions on an edge in fact has a larger symme- 
try (note that there are no interactions to consider in 
the theory). The fields are equivalent to four Majorana 
fermions, and there is S0(4) symmetry. We note that, 
as Lie algebras, S0(4)= SU(2)xSU(2), and here the first 
SU(2) can be identified with the spin-rotation symme- 



try group, while the second contains the U(l) symmetry 
just discussed as a subalgebra, generated by, say, rota- 
tion about the z axis in the second space. Thus the 
Dirac field and its conjugate can be viewed as carrying 
spin M = ±1/2 under the second SU(2). Alternatively, 
viewing the theory just as two Dirac fields, these can 
be bosonized, and we obtain two chiral bosons. The al- 
lowed "charge" states for the edge, which take values in 
the Cartan subalgebra of S0(4), lie on a two-dimensional 
lattice. This lattice is a direct product of two copies of 
the weight lattice of SU(2). Points in the lattice sim- 
ply describe the total Sz of spin and the total M on 
the edge. The same desciption applies to a vortex, since 
we can view it as an edge rolled up into a small circle. 
In the latter case, the half-flux quantum we assume in 
the vortex corresponds to changing to periodic boundary 
conditions for the fermions on the straight edge, for all 
components of the fermions. In addition to these differ- 
ent boundary condition sectors, there are also selection 
rules from the global quantum numbers. These are sim- 
ilar to the rules described in detail in Q. Specifically, 
for N even, one can have the ground state with no edge 
excitations, or one can create fermions on the edge, but 
only in even numbers. In the present situation, each of 
these fermions can be in any of the four states in the 
representation of spins (1/2,1/2) under SU(2)xSU(2). 
For odd total fermion number, there must be one un- 
paired fermion, which we can put on the edge to obtain a 
low-energy state. Then the charge (or particle number) 
differs by one from the ground state, in addition to the 
non-trivial S0(4) quantum numbers. Additional fermion 
pairs can be excited in this case also. If we consider two 
parallel edges, as for a system on a cylinder, then we 
build up the full spectrum by applying these rules to the 
two edges together, and there will be different sectors 
corresponding to the presence of either zero or one-half 
of a flux quantum threaded through the cylinder. In the 
FQHE application, there are also chiral bosons for charge 
excitations on the edge, and there are fractional charge 
sectors for each edge, though the total charge must be 
integral, and the total charge (or particle number) is cor- 
related with the fermion excitations through its parity, 
as already explained. 

The quantum numbers of the vortices (or FQHE quasi- 
particles) correspond closely to the possible quantum 
numbers for the edges, and are obtained in a similar 
way. As a simple example, we again consider two vor- 
tices. The situation is equivalent to the existence of a 
single fermion zero mode for the pair, which can be oc- 
cupied by any of four types of fermions. Thus there is a 
total of 16 states, of which half have N even, half have 
N odd. These should be analyzed as a product of four 
states possible for each vortex. The states for a vortex 
transform as (1/2, 0) © (0, 1/2) under S0(4), that is they 
either carry spin 1/2 and no M, or vice versa. These 
correspond to the states for two edges with a half flux 
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quantum through the cyHnder, in which the zero mode 
shared by the two edges can be occupied by any of the 
four types of fermions. For any even number of vortices, 
the results are similar, and the counting of degenerate 
low-energy states is fully accounted for by the four states 
for each vortex, and hence there is no nonabelian statis- 
tics. Thus a d-wave paired state, or even a superconduc- 
tor, in two dimensions, possesses vortices that may carry 
spin 1/2, but not simply because a fermionic quasipar- 
ticle can sit on the vortex. In the FQHE, these vortices 
also have well-defined charge of ±1/(2(7) for filling factor 
V — 1/q, where q must be even when the particles are 
fermions (such as electrons) . This is also true in all other 
cases discussed in this paper. It arises from the effective 
half quantum of flux that the vortices carry, and they can 
exist only in even numbers if the system has no edges. 

We hope that the above discussion gives a sufficient 
impression of the quantum numbers of the vortex states. 
For readers familiar with the general theory of abelian 
FQHE states, we will now give a precise definition of the 
structure of the state, using the Gram (or K) matrix lan- 
guage, which specifies the lattice formed by the possible 
quantum numbers (including charge) of the FQHE quasi- 
particles, as well as their statistics, and the order param- 
eters and chiral algebra of the edge theory [^;^ 78 4|] . 

The full lattice of possible quantum numbers of a vor- 
tex, or the total quantum numbers of a set of multiple 
vortices, is denoted A* (as in Ref. p8[). We will describe 
it first as a set of vectors in R"^, using an orthonormal 
basis with the standard inner product. Then A* consists 
of the set of vectors of the form 



(ri/\/2,r2/\/2,r3/[2V^) 



(49) 



where ri, r2, r^ are integers obeying ri + r2 + ^3 = 
(mod 2). Also the statistics of the excitation is O/tt — 
v^, where 6 is the phase picked up by exchanging two 
identical such quasiparticles, and the conformal weight 
of the corresponding operator in the edge theory is v^/2. 
The basis has been chosen so that the three quantum 
numbers carried by the excitations are proportional to 
the components in this basis. In fact, the spin Sz = f'i/2, 
M = r2/2, and the charge is Q = r^/{2q). From these 
rules we see that the smallest possible vortex (the one 
with the smallest nonzero v^) carries either Sz — ±1/2 
or M = ±1/2, and charge ±.l/{2q), as stated above, and 
has statistics O/tt = 1/2 + 1/(4(7) (mod 2). It is easy to 
verify that the full lattice A* is obtained as the set of 
all linear combinations, with integer-valued coefficients, 
of the vectors describing the quantum numbers of the 
basic (smallest) vortices, and thus represents the possible 
quantum numbers of any collection of vortices. 

The excitation lattice A* is the dual lattice to the 
condensate lattice A, and A is a sublattice of A*. In 
the present case, A is the set of vectors w that are lin- 
ear combinations with integer coefficients of the vectors 



(±l/-\/2, ±1/a/2, iv/9) (with all three signs independent) 
which represent the underlying particles. These represent 
possible states for a particle (electron) tunneling into an 
edge, and these electron operators, as we may term them 
in spite of the emergent M quantum number, are usually 
part of the condensate lattice, as they are in the hierar- 
chy theory. (An exception to this is the strong-pairing 
phases, where only operators of charge a multiple of two 
appear.) The fact that A and A* are dual means that for 
any v S A*, we have v • w G Z for all w S A, and vice 
versa. It suffices to check this for the w's representing 
the electron operators. 

Both lattices possess neutral sublattices, that is lat- 
tices of vectors such that Q = 0. The neutral sublat- 
tice of A, denoted A^, consists of vectors with ri and 
r2 even. Thus these form a direct sum Z Z of one- 
dimensional lattices. Each of the latter can be viewed as 
the root lattice of SU(2) in Lie algebra theory, and A-'- is 
the root lattice of SO (4). The neutral sublattice of A*, 
denoted A*-"", is the set of vectors with r^ — 0, and so 
ri + r2 =0 (mod 2). This is a sublattice of the weight 
lattice of SO (4), which would be the dual of A-*- as a 
two-dimensional lattice. The simplest nontrivial neutral 
vector is of the form (±l/\/2, ±l/-\/2, 0) (with indepen- 
dent plus and minus signs), and these represent the neu- 
tral fermions, that is the BCS quasiparticles considered 
in this paper. These cannot be created individually on 
a single edge; only excitations lying in A, such as even 
numbers of such fermions, can be p2|. 

An integral basis for a lattice is a set of vectors in the 
lattice that are linearly independent (over R), such that 
all vectors in the lattice can be written as linear combina- 
tions of those in the set, with integer coeffecients. Such 
a basis cannot be an orthogonal set of vectors, unless 
the lattice is a direct sum of one-dimensional lattices. In 
our case, a convenient integral basis for A (other than a 
suitable set of three of the electron operators above) , is 
ei = (l/V2,l/V2,Vg), 63 = (V2,0,0), eg = (0,\/2,0). 
The Gram matrix of the lattice is the set of inner prod- 
ucts of these vectors, G,;j = e^ • e^, and in this case is 



G 




(50) 



The diagonal structure of the lower-right 2x2 block re- 
flects the direct product nature of the SO (4) root lattice. 
The determinant of G, det G = Aq, determines the index 
of A as a subgroup of A* , that is the number of equiv- 
alence classes of vectors in A* modulo A. It gives the 
number of ground states of the system on a torus, or 
the number of sectors of edge states. Factoring off the 
center of mass degeneracy (7, we find that there are four, 
in agreement with the analysis based on the quasipar- 
ticle effective Hamiltonian. The Gram matrix, together 
with distinguished vectors that specify the charge Q and 
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spin Sz quantum numbers, is sufficient information from 
which to reconstruct the lattices A and A* , and hence the 
universal aspects of the phase, such as ground state de- 
generacies, quasiparticle statistics, and the theory of the 
edge states. In this context, the Gram matrix is often 
called the K matrix. 

This completes the analysis we will give of the phase. 
It is a generalized hierarchy state in the sense of Ref. 
p8[ , and resembles the 331 phase. The latter lacks the 
SU(2) of spin and hence has a two-dimensional lattice; 
its structure was described in detail in Ref. |43E3| . 

We may now compare the universal properties of this 
state with others analyzed previously. We find that sev- 
eral other constructions of this spin-singlet phase have al- 
ready been given. In Ref. ^^, this was mentioned briefly 
as the structure of both a spin-singlet state for z/ — 1/2 
proposed by Jain [g2| (see also Ref. |Q), and one pro- 
posed by Lee and Kane ||7J] . A more detailed analysis was 
given in Ref. [|73| , where it was also identified with a hi- 
erarchical construction. In the latter, one starts with the 
Halperin spin-singlet 2/5 state [Q, which is a 332 state 
in the Coulomb plasma language, and then one makes a 
finite density of spinless quasielectron excitations of that 
state, each carrying charge 2/5. The quasielectrons are 
then put in a Laughlin 1/2 state, to obtain a singlet state 
with v = 1/2. The hierarchical step implies that the re- 
sulting state has a three-dimensional lattice. In Ref. Q, 
this and the Jain construction were shown to coincide. 
Unfortunately, the formulas there contain a small mis- 
take: the final basis vector in Eq. (4.8) in that reference 
should be reduced by a factor of two, as should the en- 
tries in the first row and column of the Gram matrix in 
Eq. (4.9) there. The resulting matrix is then identical to 
G above, with q = 2, after permuting the basis vectors. 
This basis is the natural one for the hierarchy approach. 
In our G above, the top left 2x2 block (with q = 2) 
describes the Halperin 2/5 state, and reflects its origin. 
(The 2/5 state itself has the same G as the spin-polarized 
hierarchy 2/5 state, which descends from u = 1/3, as re- 
flected by the g -I- 1 = 3 at the top left.) 

The Jain v — \/q spin-singlet state was proposed as a 
trial wavefunction, namely 
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xi ^2X1,1- 



(51) 



Here Xm stands for the wavefunction for m filled Landau 
levels of spinless particles, (so xi is the Vandermonde 
determinant or Laughlin- Jastrow factor) and xu is the 
lowest Landau level filled with particles of both spins. We 
ignore the Gaussian factors in these wavefunctions, and 
have omitted the projection to the LLL. The filling fac- 
tor is again 1/g. This wavefunction can be loosely viewed 
as a Coulomb plasma of particles carrying charge, spin 
and another quantum number M = ±1/2 to represent 
the two Landau levels in X2 ■ The exponents in the wave- 
function correspond to inner products of corresponding 



vectors, which are just those of the four electron opera- 
tors with (5=1. This establishes the equivalence, as for 
the hierarchy states in Ref. Q. In fact, the extra SU(2) 
that appears here acts on the LL indices in X2 iu Jain's 
function, just like the SU(n) that appeared for n Landau 
levels in the spin-polarized Jain states p3| and the corre- 
sponding hierarchy states with n levels in the hierarchy. 
The state is a combination of the Halperin spin-singlet 
structure with the spinless composite-fermion/hierarchy 
structure. We should point out that the number of flux 
for Jain's state on the sphere is N^f, = q{N — 1) — 2, the 
same as for the HR state, which follows from our anal- 
ysis, and was noticed previously [ [72[ . This is of course 
essential in order for it to be possible to vary parameters 
smoothly to reach the transition. 

Thus the same v — \/2 (more generally, v — 1/q) 
spin-singlet phase has arisen in four different ways. We 
want to emphasize that the equivalence of the univer- 
sal long-distance properties does not mean that the trial 
wavefunctions in different approaches are the same. For 
example the paired wavefunction found here and Jain's 
above do not look alike. It may be that one is a much bet- 
ter description (has a much larger overlap with an exact 
ground state) for medium size systems than the other, 
even though they describe the same phase. The equiv- 
alence found here is analogous to that between the 331 
state and the A-phase p-wave state, however in that case 
there was an exact equivalence of certain wavefunctions 
through the Cauchy determinant identity. We may still 
expect some equivalence in the long-distance form of the 
wavefunctions. 



C. Induced Chern-Simons actions and analogs of the 
Hall conductivity 

The arguments given so far for the edge states and for 
zero modes on vortices, on which the identifications of 
the weak-pairing phases have been based, may appear 
not to be very well-founded, as they have been based on 
analyzing the BdG equations for special forms of the gap 
function and variation of the parameter fi, though we did 
argue by continuity that the states found must persist as 
the equations are varied while staying in the same phase. 
In this Subsection, we argue that the results we have ob- 
tained are in fact very robust, because the edge states, 
and the form of the bulk theory described by the Gram 
matrix of the condensate lattice, are required as a conse- 
quence of transport properties of the bulk weak-pairing 
phases. These transport properties are the quantization 
of the spin and heat analogs of the Hall conductivity, 
which we prove explicitly for the spin case. They imply 
the existence of certain edge state structures, just as in 
the case of charge transport in the QHE, and when the 
weak-pairing phases correspond to abelian FQHE states, 
the Hall spin conductivity is actually a part of the Gram 
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matrix description. The remainder of this Subsection 
discusses these points, but the technical details are rel- 
egated to the Appendix. We should point out that the 
form of the argument has already appeared in Ref. | |57| ] 
for the cases with SU(2) or U(l) spin symmetries, though 
the explicit derivation of the Hall conductivities was not 
given there. 

For the cases of spin-singlet pairing, and of p-wave 
pairing with an unbroken U(l) symmetry, we derive in 
the Appendix the Hall spin conductivity a^y, and show 
that in any fully-gapped translationally-invariant super- 
conducting phase it is given by a topological invariant, 
which within the BCS approximation is proportional to 
the same winding number we discussed earlier in Sec. 
U and subsequently. Some similar statements appeared 
earlier, but were for the charge Hall conductivity [§,||, 
and also for the existence of a Hopf term and of a non- 
abelian CS term, in the two-dimensional A-phase with 
SU(2) symmetry broken spontaneously |5[|. In this paper 
we restrict ourselves to the Hall conductivies for con- 
served quantities for which the corresponding symmetry 
is not broken spontaneously in the paired state, such as 
spin and energy. For quantities for which this is not 
true, such as charge, we do not find quantized Hall con- 
ductivities in superconductors (though of course we do 
find quantized Hall conductivity in the QHE systems). 
Since our point of view differs somewhat from that in, 
for example, Ref. B, we give a self-contained discussion. 

To be precise, for the Hall spin conductivity in a spin- 
singlet paired state, we find that if the particles are 
viewed as carrying spin 1/2 (we usually set h = 1), the 
Hall response to an applied spin analog of the electric 
field, such as a gradient in the Zeeman splitting, is 
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where m is the winding number, which is ±2 in a d-wave 
weak-pairing phase. We have written the Hall spin con- 
ductivity in this form to emphasize the similarity to the 
usual e'^/h, with e replaced by h/2 here. In h = 1 units, 
we obtain in our d-wave weak-pairing phase a^y = l/47r. 
We chose the -I- sign, since in the FQHE applications of 
the / — —2 state, the edge modes propagate in the same 
direction as the charge modes. This agrees with Refs. 
[pq,M, where different arguments were used. For the p- 
wave case, we view the fermions as carrying isospin ±1, 
and hence the Hall spin conductivity we obtain is 
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where the winding number to is ±1 in the weak-pairing 
p-wave phases, again with m = I for the ^ = — 1 case. 

These results agree with the descriptions we have al- 
ready given of the weak-pairing FQHE phases using the 
Gram matrix or lattice theory. In fact we should point 



out that a quantized Hall spin conductivity is not unusual 
in FQHE systems, though it is not always emphasized. 
It occurs for example in any spin polarized state, such 
as the integer and Laughlin states with ly = 1/q, because 
the electrons carry spin 1/2 as well as charge, and so the 
two Hall conductivities are proportional. It also occurs 
in some spin-singlet Hall states (abelian examples were 
discussed in [^), including the v = 2 state with the 
LLL filled with both spins, and the Halperin mmn states 
with m = n + 1, m odd, which are a generalization of 
the 1/ = 2 state. In these cases, we obtain the full SU(2) 
version of the Hall spin conductivity, with a^y taking the 
same value as in the d-wave weak-pairing phase. The 
same quantized Hall spin conductivity was also found in 
certain spin-liquid states for lattice antiferromagnets [£9| . 
The wavefunctions of these states are the same as that of 
the Halperin state, with the charge degree of freedom re- 
moved (i.e. the wavefunction is 1/2, 1/2, —1/2), and the 
spin-1/2 particles restricted to lattice sites, explaining 
this result. 

The result that the Hall conductivity, in units of (quan- 
tum number)'^//i, is a topological invariant given by an 
integral over k space, is similar to one form of the or- 
dinary charge Hall conductivity, found originally for a 
noninteracting periodic system with a rational number of 
flux quanta per unit cell, as an integral over the Brillouin 
zone [ pO[ . In systems where k is not a good quantum 
number, such as the same system with irrational flux, or 
when disorder is present, or in other geometries, includ- 
ing those with edges, that break translational invariance, 
such a topological invariant is apparently not available. 
Yet the idea of quantization as resulting from the conduc- 
tivity being a topological invariant that measures an in- 
trinsic local property of the ground state seems too good 
to give up. For the noninteracting QHE with disorder, 
the topological invariant has been extended using non- 
commutative geometry, so that the Hall conductivity is 
an integral over a "noncommutative Brillouin zone" , and 
in this way quantization has been proved even for the 
physically relevant case of a nonzero density of localized 
states at the Fermi energy ^l[. It would be interesting 
to extend this to other cases, including the paired states 
with disorder, which we discuss in the next Section. In 
the Appendix, we give arguments that the form we ob- 
tain is exact to all orders in the interactions, but only for 
a translationally-invariant system. 

Now we can use the results on the Hall spin conduc- 
tivity to argue that the edge state properties we have 
obtained arc indeed correct. In the Appendix, we derive 
the Hall spin conductivities by obtaining the induced ac- 
tion for an external gauge field that couples to the spin 
or isospin. The actions that result (cf. Ref. p3]) are 
CS terms for an SU(2) gauge field in the spin-singlet (d- 
wave) case, and for a U(l) gauge field in the p-wave case. 
Now, using either the Hall conductivities and arguing as 
in the QHE literature p3] (see also p4,p3|), or using the 
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induced actions and arguing as in the field theory liter- 
ature iQ (quite similar arguments appear in Ref. |8^] 
and references therein), we can conclude that on a do- 
main wall between phases with different cr^^'s (one of 
the phases might be the vacuum outside an edge, with 
a^ — 0) there must be chiral edge excitations. In the 
presence of a uniform spin-electric field, a spin-current 
is induced in a region with a^y ^ 0, and the normal 
component of this at a domain wall has a discontinu- 
ity, representing a net inflow of spin onto the wall. To 
avoid violating the continuity equation for the spin den- 
sity and spin current density, there must be chiral modes 
on the wall, and a "gauge anomaly" in the conservation 
of spin on the wall alone. The tangential field induces 
a nonzero divergence (i.e., an anomaly) of spin current 
along the wall, which cancels the net infiow from the 
bulk. Such an anomaly can occur only if gapless chiral 
excitations exist on the domain wall. The minimal chi- 
ral theory required to produce the anomaly is the usual 
chiral Luttinger liquid (or chiral Gaussian model) in the 
U(l) case, with the value of the couphng that corresponds 
to a free chiral Dirac (Weyl) fermion in 1 -I- 1 dimensions, 
and the SU(2) chiral Wess-Zumino-Witten (WZW) the- 
ory |8|] with k ^ m/2 = 1 in the SU(2) case. All of 
this applies even within the BCS mean field framework 
we used before, and then the edge excitations must be 
free fermions. In the p-wave case, we therefore expect 
simply a single chiral Dirac fermion to propagate on the 
edge. In the d-wave case we must have an SU(2) doublet 
of chiral Dirac fermions, which can be represented by the 
k — 1 chiral WZW model, together with an additional 
chiral U(l) degree of freedom, which we argued earlier 
must exist, and called M. We have therefore reproduced 
the claimed results about the edge states, neglecting the 
charge degree of freedom. 

We note that, when formulating such arguments for the 
nature of the chiral edge theories, we can presumably as- 
sume that the theories are unitary, conformal fields with 
local current operators for physical conserved quantities 
(as usual, Lorentz invariance may be spoiled by the pres- 
ence of different velocities for different excitations, but 
this will not matter for the statistics and other univer- 
sal properties in which we are interested). Previously, 
we might not have assumed this, because of the example 
of the Haldane-Rezayi (HR) state in particular. But we 
have learned that the HR state is at a critical point, and 
previous discussions of the edge and quasiparticle prop- 
erties of that state are irrelevant. Thus, with the demise 
of the HR state as a bulk phase, it becomes attractive to 
believe that the edge theories of incompressible FQHE 
phases are always unitary conformal field theories. With 
this assumption, in theories with SU(2) Hall spin conduc- 
tivities, unitarity of the edge theory requires quite gen- 
erally that k be an integer [p8 39|, and so a^ — k/{4n) 
in the above conventions. 

The preceding arguments do not apply to the spinless 



p-wave case, in which there is no continuous unbroken 
symmetry. This is unfortunate in view of the great in- 
terest in the nonabelian properties of the weak-pairing 
phase. But there is another Hall-type conductivity, which 
exists in all cases, including those without a continuous 
symmetry. This is the Leduc-Righi (LR) conductivity, 
which is the xy component Hxy of the thermal conduc- 
tivity, and is of course related to the transport of en- 
ergy, a conserved quantity. Like the Hall conductivities 
for charge and spin, this is a non-dissipative transport 
quantity that can only be nonzero when parity and time- 
reversal symmetries are broken. In systems with a gap 
for all bulk excitations and with chiral edge excitations, 
it can be argued that the LR conductivity is nonzero 
^ , p7| . Thus this applies to QHE systems, and to su- 
perconductors (paired systems) if there is no gapless col- 
lective charge mode. The value of Kxy can be obtained 
^ , p7| by considering a sample with two edges and a 
small temperature difference between the edges (strictly 
speaking, the following argument yields the LR conduc- 
tance, not the conductivity). The chiral excitations on 
each edge are excited by different temperatures, and this 
produces a larger heat current on one edge in one direc- 
tion than that on the other edge in the other, and hence 
a net heat current. This shows that the current is related 
to the heat capacity of each mode on an edge, times the 
velocity of the mode, summed over modes. It is known 
that the heat capacity for each mode is related to the 
Virasoro central charge c in the conformal field theory of 
the edge excitations [|o|. In n^y, the velocities cancel, 
and the remaining number is proportional to the total 
central charge of the edge theory (including the charge 
modes in the FQHE); precisely. 
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This is the fundamental relationship governing the LR 
conductivity in all quantum Hall problems. Here we as- 
sumed that all modes on an edge propagate in the same 
direction; if not, then c should be replaced by the differ- 
ence of the central charges for the right and left-moving 
theories. Related to this, the two-probe thermal conduc- 
tance of such a system is also equal to K^y (for a case with 
all edge modes propagating in the same direction), just 
as the two-probe conductance is equal to the quantized 
(Txy of the bulk. It has been found that the two-probe 
thermal conductances of the Laughlin states at various 
ly = 1/q are independent of q, even though the theories of 
the chiral Luttinger liquids contain q as a parameter [ pi| . 
This is because the central charge is c = 1, independent 
of q. 

The formula for K^y is similar in structure to that for 
the Hall spin conductivity. However, the central charge 
c does not have to be an integer, and indeed for a single 
Majorana fermion, c = 1/2 |92|, so that for free fermions 
in general, c is a multiple of 1/2. Nonetheless, we do 
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expect it to be quantized, in the sense of invariant under 
small deformations of the theory (including adding weak 
disorder). 

The argument above for the LR conductivity made use 
of the edge states. However, we want to use it to back 
up our results on the fermionic edge modes of the MR 
state, in analogy with the arguments for the other cases, 
which used the Hall spin conductivities. Hence we need 
an independent argument for the value of the LR con- 
ductivity. We believe that it should be possible to derive 
such a result, analogously to the Hall spin conductivities, 
by considering the system in external gravitational fields. 
Here the Christoffel symbols, or the spin connection, play 
the role of the external gauge fields we used in the spin 
case, but should be viewed as determined by the metric 
of spacetime, which we treat as the independent variable 
and set to the usual Minkowski metric after calculating 
responses. The role of the currents, to which the vec- 
tor potentials couple, is played by the energy-momentum 
tensor, which includes the energy flux among its compo- 
nents. The significance of changing the metric should be 
clear if we recall that equilibrium systems can be repre- 
sented in imaginary time, with the imaginary time direc- 
tion periodic, the period being l/fc^T. Thus a tempera- 
ture gradient might be viewed as changing the geometry 
of space and (imaginary) time. A more rigorous deriva- 
tion would avoid imaginary time, but should still involve 
the response to changes in the metric. The leading term 
in the induced action for the external gravitational field 
in 2 + 1 dimensions will in general be the gravitational 
CS term, which can be written in terms of the spin con- 
nection in close analogy to the usual CS terms: 

1 c /■ 2 
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where we view the spin connection w as a 3 x 3 matrix, 
determined by the metric. (Such a term was also pro- 
posed earlier for He'^-A films in Ref. [Q.) The coefficient 
contains c/24, which shows the relation to the central 
charge c in a corresponding chiral conformal field theory 
on a boundary |3Sf| , is needed to cancel the anomaly in 
energy-momentum conservation on the boundary by an 
inflow (LR "Hall" current) from the bulk |g^, as for the 
charge and spin Hall conductivities. 

We have not completed a calculation (analogous to 
those in the Appendix) of the induced action or LR con- 
ductivity for the bulk from first principles. It would 
involve coupling the underlying system to an arbitrary 
metric, then using Ward identities to relate the response 
to a topological invariant. However, if we assume a spec- 
trum of relativistic fcrmions at low energies, with a min- 
imal coupling to the gravitational field, then the calcu- 
lation can be carried out, and is known in the literature 
p3[ . The coefficient of the induced gravitational CS term 
again has the form of the same topological invariant as in 
the Appendix, but integrated over only half the sphere. 



as in the similar treatment of the spin cases. For a single 
Majorana, this yields the above form with c replaced by 
±1/4 on the two sides of the transition. Inclusion of the 
contribution of a Pauli-Villars regulator then produces 
c = on one side, c = 1/2 on the other. We expect 
that for our nonrelativistic system, where the large k be- 
havior is explicitly known, we would obtain such a result 
also, with c in a general system of paired fermions being 
proportional to the same topological invariant (winding 
number) we have seen already, and (allowing for factors 
of two associated with spin degeneracy) we would find 
c = 1/2 in the spinless p-wave weak-pairing phase, and 
also c = 1 (2) in the triplet p-wave (singlet d-wave) cases. 
(These are the results for the paired fermions, and in the 
FQHE would have to be supplemented by adding c — 1 
for the charge degrees of freedom.) Therefore we believe 
that the existence of the edge states and vortex degen- 
eracies in the MR phase can be placed on a firm footing. 
A somewhat related issue is to obtain effective ac- 
tions describing the weak-pairing phases. We emphasize 
that the CS actions discussed in this Subsection are in- 
duced actions for external fields which act as sources, and 
should not be confused with effective actions, which con- 
tain fields that should be functionally-integrated over, 
and represent the dynamics of the system at low ener- 
gies and long wavelengths. For the abelian FQHE states, 
the effective actions fall into the framework of the known 
theory, based on the Gram (K) matrix |£8[ . For the MR, 
and other similar, phases, something different is required. 
For the abelian phases, there is an evident similarity be- 
tween the induced and effective actions. For example, 
in the d-wave case, neglecting the charge degree of free- 
dom, one would expect the bulk effective action to be an 
SU(2)xU(l) CS theory, with fc = 1 for the SU(2) part 
(here, in the effective action, k must be an integer to 
maintain gauge invariance). This theory is determined 
by the requirement that it produce the correct edge the- 
ory Q. (The U(l) part could possibly be extended to 
a second SU(2) with A; = 1, to agree with the S0(4) 
edge theory discussed above.) It is known that such an 
effective theory also produces the desired induced grav- 
itational CS term in all cases |g^. By analogy, we are 
led to expect that the MR phase, where the edge theory 
involves only the c = 1/2 representations of the Virasoro 
algebra, can be described by an effective theory which 
is some sort of gravitational CS theory (similarly, there 
were earlier proposals for gravitational-type effective ac- 
tions for three-dimensional paired states M). Quanti- 
zation of such a theory should yield Virasoro conformal 
blocks in the same way that quantization of CS theory 
yielded current algebra blocks [p9|, and thus be closely 
related to the wavefunctions discussed in Sec. O and in 
|l4,E4|. This hope is encouraged by the identification 
of the parameter A and the vector potential A, which 
should be functionally integrated over in the full treat- 
ment, as the vielbein and spin connection of 2-t- 1 gravity. 
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or at least as the part relating to S0(2) rotations of space 
only, as discussed in Sec. |l|. The hope of producing con- 
formal blocks in this manner from a gravitational analog 
of CS theory has been around for a long time |39|, but 
does not seem to have reached fruition, in spite of an in- 
teresting attempt by Verlinde BJ]. Such a theory would 
be an interesting, possibly more natural, alternative to 
the "conventional" approach, along lines anticipated in 
Ref. Q , of a CS analog of a coset construction Q . 

We also wish to comment on whether our results imply 
that fractional and nonabelian statistics occur in paired 
superfluids, as opposed to FQHE systems. For example, 
the weak-pairing d-wave and spin-triplet p-wave phases 
have nontrivial Hall spin conductivity, and the smallest 
possible vortices carry spin (or M) quantum numbers. 
Hence the spin degrees of freedom contribute a fractional 
amount to the Berry phase on exchanging such vortices. 
Similarly, exchanging vortices in the weak-pairing spin- 
less p-wave phase should produce a matrix action on the 
space of degenerate states we have identified, which we 
may be tempted to term nonabelian statistics. However, 
although these contributions from the spin (or fermion 
number) sector are well-defined, in a neutral superfluid 
the charge degree of freedom is gapless, and the vortices 
act on the charge (particle number) variables also, so 
as with vortices in a simple neutral superfluid, which 
carry no well-defined particle number, the contribution 
to the total Berry phase is not well-defined, due to the 
charge sector [B5| . Nonabelian statistics is still meaning- 
ful, modulo phase factors. In the incompressible FQHE 
phases, this problem disappears, and the statistics prop- 
erties have been characterized above, in detail for the 
abelian cases. Also, in a superfluid with a Coulomb in- 
teraction, there is again no problem, even if the interac- 
tion is '^ 1/r, which does not produce a plasmon gap in 
two dimensions. Vortices are neutral because of screen- 
ing, and so nonabelian statistics, or fractional statistics in 
the triplet p-wave case, can occur, with no contribution 
from the charge sector. There are also "neutral" vortex 
excitations with no net (spin-independent) vorticity act- 
ing in the charge sector, which would not be subject to 
the problem in compressible superfluids, but these are 
found not to have fractional statistics. Note that a simi- 
lar problem as for the charge sector (in the compressible 
case) occurs in the spin sector in connection with any 
subgroup of SU(2) that is spontaneously broken. 



V. EFFECTS OF DISORDER ON THE 
TRANSITIONS 

In this Section, we discuss the effects of disorder on 
the phases above, and on the transitions between them. 
We consider the phase-coherent, zero-temperature case, 
and neglect all interactions between the quasiparticles 
(including the gauge field fluctuations). The problem 



then reduces again to the quasiparticle effective Hamil- 
tonians Kgg, this time with ^k and Ak replaced by op- 
erators that are local in position space, with short-range 
correlations of the disorder. We consider the problems 
above in reverse order, starting with the d-wave case, 
which has the most symmetry (SU(2) of spin rotations), 
then d- or p-wave with only U(l) symmetry, and finally 
the spinless p-wave case, with no continuous symmetry. 
The first case has been studied recently in the context of 
disordered superconductors |p3-pq], while we will argue 
that the second maps onto the usual noninteracting QHE 
transition, and the last includes an unusual intermediate 
phase where the heat transport is similar to that in a 
disordered metal. In all cases, we expect the qualitative 
results to be unaffected by interactions (or quantum fluc- 
tuations around the mean field theories used) , though the 
universality classes may be changed. 

The problem of the noninteracting BCS quasiparticles 
in disordered paired fermion systems was discussed by Al- 
tland and Zirnbauer (AZ) [^, where a symmetry classifi- 
cation of random matrix ensembles was proposed, that is 
analogous to the familiar classification for ordinary one- 
particle Hamiltonians into orthogonal, unitary, and sym- 
plectic ensembles (or symmetry classes). In these ensem- 
bles, no particular value of the energy is singled out as 
special, so that they apply to phenomena near generic 
energy (or Fermi energy) values, where the average den- 
sity of states is nonzero. The disordered paired systems 
differ first in that the fermion number is not conserved, 
because of the pairing terms. The full Hamiltonian of 
course conserves number, but this involves the collective 
response of the condensate. In applying these models in 
the FQHE, this is again true, but leads to the Hall re- 
sponse, not a superfiuid response. Therefore we disregard 
particle number transport, and concentrate on conserved 
quantities carried by the quasiparticles only. The latter 
quantities include spin, when this symmetry is not bro- 
ken spontaneously by the pairing. Then the classification 
of ensembles by AZ is according to whether time-reversal 
symmetry (T) is broken or not, and whether or not there 
is an unbroken SU(2) spin rotation symmetry; the lat- 
ter is unbroken in spin-singlet paired states. By making 
certain transformations, such as a particle-hole transfor- 
mation on the I spins in the cases where Sz (at least) 
is conserved, the quasiparticle Hamiltonians can be re- 
lated to number-conserving Hamiltonians, and thus to 
single-particle or random matrix problems. In this way, 
AZ identified four classes of random-matrix ensembles 
for disordered paired systems, which they labelled C, CI, 
D, DHL For the quasiparticle Hamiltonians, the zero of 
quasiparticle energy E is a special point, where in most 
cases the average density of states of the quasiparticles in 
the disordered system vanishes. Thus these four classes 
are distinct from the usual three mentioned above. 

The case of spin-singlet paired states with disorder 
was considered in more detail ||5l|,p3|,p4[ . The symme- 
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(a) class C 
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(b) class D 



FIG. 2. Proposed renormalization group flow diagrams for 
(a) the unitary ensemble (IQHE), as in pa|, and class C, 
and (b) class D. The dashed curves represent schematically 
the (nonuniversal) bare values of the coupling parameters. 
Other features are universal when the renormalized couplings 
are deflned using the as-measured conductivity parameters, as 
explained in the text, and repeat periodically in the 9 variable. 



try classes are CI (with time reversal symmetry), and C 
(without time reversal). In the particular case of class C, 
it is natural to consider in two dimensions the possibility 
of a Hall spin conductivity crj„, as in the weak-pairing 
d-wave phase [p6| , p7[ . There are nonlinear sigma model 
formulations for class C, either using replicas, which lead 
to target manifold Sp(2n)/U(n) (with n — > 0) in the com- 
pact formulation p3[ , or supersymmetry, which leads to 
the target supcrmanifold Osp(2n|2n)/U(n|n) (with n> Q 
arbitrary) [ p2[ . The Hall spin conductivity shows up in 
that a topological term with cocfhcient proportional to 
at,, can be included in these models, when the dimension 
of space is two |53] . Analysis of the model, and numer- 
ical work on a network model with the symmetries ap- 
propriate to class C pq] , has shown that this system has 
phases with quantized aj which are the same, or multi- 
ples of, those we discussed in the pure system in Subsec. 



[VC. These phases can therefore be viewed as the disor- 



dered analogs of the weak- and strong-pairing phases. In 
these phases, the BCS quasiparticle (fermion) excitations 
at low energies have a non-zero density of states which 
vanishes quadratically as -E -^^ 0, and these states are 
localized [Q. The situation is thus similar to the usual 
QHE, in which the low-energy excitations are localized 
and have nonvanishing density of states at the Fermi en- 
ergy. Localization is necessary to obtain the quantized 
Hall conductivity. In the weak-pairing phase, we there- 
fore expect that the results for the edge and statistics 
properties obtained in Sec. I\ are still valid when disor- 
der is included. We do not expect this conclusion to be 
affected by the inclusion of interactions in the analysis. 

The transition between these two phases, and the role 
of a^y, can be understood via a renormalization group 



(RG) flow diagram, shown schematically in Fig. g (a). 
The flows can be thought of as representing the values of 
the local conductivity parameters a%^, a^y that would be 
measured at a given length scale, and how they change 
with this scale. In the nonlinear sigma model, similar 
parameters appear in the action of the quantum field 
theory. It is natural to define the renormalized values 
of these couplings at any scale to be the conductivities 
that would be measured at that scale, in which case the 
flows are the same as the RG flows within the field the- 
ory model. The conductivities should be understood in 
this way from here on, instead of as the bare values, with 
which they coincide only when the scale is of order the 
mean free path. Quantized values refer to the renormal- 
ized values at very large scales. In the case of class C, the 
form of the flows is identical to that in the IQHE for non- 
interacting electrons |96|. Similarly to the usual IQHE, 
the transition occurs because of a nontrivial fixed point, 
at which ct* is midway between adjacent quantized val- 
ues, and CT^j. will take some nontrivial universal value. 
However, the spin quantum Hall transition in class C is in 
a different universality class from the IQHE. Recent nu- 
merical work has obtained some of the critical exponents 
for this transition, both from the network model p6[ , and 
from a mapping of the network model to a supcrsymmet- 
ric vertex model and a superspin chain [p5 58 5^1 , which 
was then analyzed numerically p7| . The results are in ex- 
cellent agreement with exact values of some exponents, 
which were proposed using the relation of the supersym- 
metric vertex model to classical percolation psi. For 
example, the localization length exponent in v — 4/3, 
and at the critical point the density of states vanishes 
as £'^/^. The effect of interactions on these results is 
presently unknown. 

Inclusion of a Zeeman splitting ex h, which was ne- 
glected so far, will split the transition, and the phase 
diagram as a function of ^ and h will be similar to Fig. 4 
in Ref. [p7| . There will be an intermediate phase in which 
ctJ„ is quantized and halfway between the quantized val- 
ues in the phases on either side. (This is somewhat like 
the intermediate MR phase produced by tunneling t in 
the p-wave case.) In the present case, the Zeeman term 
leaves unbroken a U(l) subgroup of the SU(2) symmetry 
present in the spin-singlet paired state; for h parallel to 
the z axis in spin space, this U(l) is generated by Sz. The 
pairing is still between spin | and spin |, and so the sym- 
metry is the same as in the p-wave case: with disorder, 
the distinction between p-wave and d-wave, and between 
spin-singlet and spin-triplet, is lost. Therefore we expect 
that in this intermediate phase, there is a single chiral 
Dirac fermion mode on the edge, and that the statistics 
properties are the same as those of the 331 state. Also, 
the transitions from the weak-pairing abelian phase, in 
which the statistics properties will be unchanged even 
though the SU(2) symmetry is broken (similar to the ef- 
fect of t on the 331 state), to this phase, and from this 
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phase to the strong-pairing phase with a. 
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pected to be in the universahty class of the usual IQHE 
|p6|"p8[. Here the unbroken U(l) of spin is playing the 
role of the particle number in the usual IQHE case; the 
real particle number is of course still not conserved in the 
paired state. The appearance of this symmetry class in 
a disordered superconductor was seemingly overlooked 
by AZ; in fact, there are two such possible classes, in 
which there is an unbroken U(l) [not SU(2)], and T may 
be either unbroken or broken. By applying the methods 
used by AZ in the case of broken T and unbroken U(l) 
symmetry, one is led back to the usual unitary ensemble, 
with the spin "f and spin J, quasiparticles playing the role 
of particles and holes, respectively. It follows that these 
phases at non-zero h possess a nonvanishing density of 
(localized) states for the quasiparticles at -E ^ 0. In the 
two-dimensional case, the unitary class admits a topo- 
logical term with coefficient oc cr^ in our case, and thus 
an IQHE transition. This is consistent with the results 
in Rcf. |5g]. Note that the relevant interaction effects in 
the case of paired states and the IQHE may turn out to 
be different, however, so the equivalence might not hold 
when interactions are taken into account. 

Turning to the p-wave states relevant to the double- 
layer system in the FQHE, for t = and with non- 
zero disorder we are once again in the situation just dis- 
cussed of broken T and unbroken U(l). Hence for this 
case, we again expect, within our model of noninteracting 
quasiparticles, a transition between quantized a^y phases 
(analogs of the weak and strong-pairing phases) that is in 
the universality class of the noninteracting IQHE. With 
t 7^ 0, we saw in the pure case that the transition splits 
into two, each in the class of the spinless p-wave case, 
and the intermediate phase had the properties of the MR 
state. The non-zero tunneling breaks the U(l) symme- 
try. We consider this case next, and then return to its 
application to the double-layer system. 

The symmetry classes of pairing of spinless fermions, 
or with SU(2) symmetry fully broken by the Hamilto- 
nian, are denoted Dili (with unbroken T) and D (with 
broken T) by AZ. These cases were not analyzed in two 
dimensions previously, but some of our results have been 
found independently in Ref. |97[ (see also a remark in 
Rcf. iQ). Our interest is in class D, with broken T. 
In this case the nonlinear sigma model target (super- 
) manifold is S0(2n)/U(n) (using replicas in the com- 
pact formulation, with n -^ 0), or Osp(2n|2rt)/U(n|n) 
using supersymmetry (with n > arbitrary) p^]. In the 
supersymmetric formulations, classes C and D differ in 
that, while the bosonic submanifolds are real forms of 
SO(2n)/U(n)xSp(2ri)/U(7i) in both cases, in class C the 
first factor is noncompact and the second compact, and 
in class D it is the other way round (corresponding to the 
compact replica formulation) |52| . (These statements are 
for n > 1. For n = 1, the first factor is a single point in 
both cases.) In the two-dimensional case, class D admits 



a topological term in the nonlinear sigma model, like class 
C. In the case of class D, there is no continuous symmetry 
in the underlying fermion problem. The only candidates 
for the physical meaning of the couplings in the nonlin- 
ear sigma model are in terms of thermal conductivities, 
since energy is still a conserved quantity p3| . The diago- 



nal (dissipative) thermal conductivity Kx 



•^yy, 



and the 



off-diagonal LR conductivity n^y have the dimensions of 
k^T/h {h is Planck's constant, = 27r in our units) at low 
temperatures. We define 






(56) 



where the numerical factor of 7r^/3 is th at wh ich arises 
in the quantized values of n^y (see Sec. IV C ), and so 
may be conveniently included here. Then we expect that 
the sigma model couplings are (similarly to the charge 
and spin transport cases, where no factors of k^T are in- 
volved) l/g^ ^ iixxi and 9 — Airk^y. Since the quantized 
values of hr^y are k^y — c, which is a multiple of 1/2 for 
Majorana fermions, we have arranged that the quantized 
values of 6 would be multiples of 27r. We expect that if 
the nonlinear sigma model for class D is derived micro- 
scopically at weak coupling, which wc will not attempt 
here, then the above relation for 9 will hold. Wc note 
here that in class D, the density of states at i? ^ in 
the localized (quantized k^y) phases is expected to ap- 
proach a nonzero nonuniversal constant, as one can see 
from the random matrix expressions in AZ p2| , using 
an argument in Ref. Ba]. We want to emphasise that 
the nonlinear sigma model in class D describes only the 
case of very generic disorder, and not necessarily more 
restricted forms of disorder. We will return to this point 
below. 

We may now consider the form of the RG flows for class 
D in two dimensions. We begin with perturbation theory 
at weak coupling, that is large k^x- This was considered 
previously for classes C and CI 1 51 5^ , and for class D is 
mentioned in Ref. ||5J]. We define the RG beta function 
as 



Pxxig 



rflnL' 



(57) 



where g^ is the nonlinear sigma model coupling squared, 
and g^ ~ l/k^x, and L is the length scale on which the 
renormalized coupling is defined. In two dimensions this 
has the form 



Pxxig') = ag^ + 0{g') 



(58) 



at weak coupling (small g^); effects of the topological 
term involving 9, if the model admits one, are nonper- 
turbative in 5^, of order e"**/^ , and contain the 9 de- 
pendence. (Also, there will be another beta function 
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Pxy = d9/d\nL for 9, which will be entirely nonper- 
turbative.) Here a and b are constants. In classes C 
and CI, a is positive |5l|,|5^, and the flow is towards 
strong coupling, that is localization, as shown in Fig. 
(a) for class C. However, in classes D and Dili, a 
is negative, and in fact equal to minus its values in C 
and CI respectively |^,Q. The reason hes in the re- 
lation of the target (super-)nianifolds in the nonlinear 
sigma models, described above. This relation is similar 
to that between the manifolds in the symplectic (spin- 
orbit scattering) and orthogonal (potential scattering) 
ensembles of the usual random matrix or localization 
problems, where the a's also have opposite signs. The 
origin of this is that the a term in the beta function of 
any nonlinear sigma model is related to the Ricci curva- 
ture tensor of the target manifold. When we compare 
these for the compact and noncompact versions of the 
"same symmetry", such as Sp(2ri)/U(n) (at n a positive 
integer), we find that they are of opposite sign: the non- 
compact case has negative curvature, the compact pos- 
itive. These geometric phenomena for symmetric Rie- 
mannian (non-super-) manifolds are discussed by Hel- 
gason Q. The noncompact space Sp(2n, R)/U(n) at 
n ^ represents class D in perturbation theory, as does 
the compact space S0(2n)/U(n) in the same limit, and 
so their perturbative /? functions are equal. This estab- 
lishes the result using replicas. Likewise for the target su- 
permanifolds, the factors in the bosonic parts each have 
opposite curvature in the two cases, and this presum- 
ably extends to the full supermanifold. Consequently, 
a has the opposite sign in the two cases. A similar re- 
sult also holds in the principal chiral models with target 
spaces Sp(2n) and S0(2n) in the compact replica ap- 
proach, with n ^ 0, which describe classes CI and Dili 
respectively. Note that in the unitary case the target 
supermanifold is \J{n,n\2n)/[U{n\n)xlJ{n\n)], and the 
bosonic part contains just the compact and noncompact 
forms of U(2n)/[U(n)xU(7i)], so the model maps to it- 
self under interchange of compact and noncompact, and 
hence the net coefficient a in the unitary case vanishes, 
as is well-known. 

We now try to find the simplest possible flow diagram 
compatible with the weak-coupling behavior and some 
other simple requirements. In the unitary case, the flows 
in Fig. |2| (a) can be considered to be the simplest pos- 
sible, on including nonperturbative 0-dependent effects 
that cause the attractive flxed points to be at = (mod 
27r). These flows are in fact obtained if one takes the 
dilute instanton gas calculations of Pruisken, which are 
the leading nonperturbative effects at weak coupling, to- 
gether with the perturbation theory result just discussed, 
and uses these forms for all coupling values |9q|. The 
nontrivial fixed point at 6 = tt and at some universal 
axx controls the transition between the quantized fixed 
points in this case. The picture obtained from these flows 
seems to be in qualitative agreement with what is known 



from numerical work for this transition in the unitary 



ensemble, and for class C |56 



For class D we can try to guess the nonperturbative 
form of the flows without calculation. In view of the 
weak-coupling result, we could try reversing the arrows 
on the flows for unitary and class C. However, we also 
expect that the stable, attractive fixed points, which will 
represent quantized values of kxy at kxx = 0, will be at 
6 = (mod 27r) again, not at tt (mod 27r). In particular, 
this means that an insulating phase with quantized kxy = 
is possible. It seems reasonable that sufficiently small 
Kxx can produce localization when kxy = 0, in spite of 
the fiow to kxx = oo in the weak-coupling region, just as 
in other localization problems, including the case of spin- 
orbit scattering, and this should be stable against small 
changes in kxy In order to achieve this, we also shift 
the flows by it along the 9 direction. The result is the 
flow diagram shown in Fig. y (b). The interesting non- 
trivial fixed points now occur on the lines 9 = (mod 
2tt). These flows could be checked in the weak-coupling 
region by comparing them with a dilute instanton gas 
calculation as in the unitary case. Indeed, if the latter 
calculation is assumed to give the same form as in the 
unitary case pq] , as is plausible, then the competition 
with the perturbative terms will give the flows as shown. 

In order to use the RG flow diagram to make predic- 
tions about the effects of disorder on the pure transi- 
tion from weak to strong pairing with T broken and no 
spin-rotation symmetry, we need to know where the bare 
values of kxx and kxy lie on the diagram. In the usual 
IQHE unitary case, and also for class C, the values are 
shown as the dashed curve in Fig. pi (a). If one uses the 
self-consistent Born approximation to obtain the values 
in the IQHE case, for disorder weak compared with the 
cyclotron energy, then one obtains a semicircle in the axx- 
(Jxy plane [|96| . The precise position of the curve is unim- 
portant, but it associates the transition, at which Uxy 
is half an odd integer, with the middle of the disorder- 
broadened Landau bands. Similar behavior occurs for 
class C. 

In the present case of class D, we again expect the bare 
values of kxx f^nd kxy to lie on an arc, as shown in Fig. 
(b). These values pass through the quantized points 
at kxx = {g^ = oo) and = (mod 27r) (quantized 
kxy)- This is reasonable, as these are the values in the 
two phases in the pure case, and disorder that is weak 
compared with the gap |/i| in the spectrum should pro- 
duce only small corrections to these values. Connecting 
these regions with the dashed arc, we always produce a 
curve of the form shown for topological reasons. This 
curve intersects the separatrices shown, which fiow into 
the non-trivial fixed points. We see that the regions near 
the quantized fixed points fiow into those fixed points, so 
that the quantized phases, one of which corresponds to 
the nonahelian statistics phase in the pure case, still ex- 
ist in the presence of disorder according to our proposed 
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flow diagram. But there is an intermediate set of bare 
values near = tt (mod 27r), which flow to weak coupling, 
and at large scales they map onto the entire interval of 
kxy values between the two quantized values in question, 
with a kxx that increases logarithmically with L, accord- 
ing to the weak coupling beta function /3xx above. This 
is therefore an intermediate phase with metallic behavior 
of the thermal conductivities, between the two quantized 
phases. The intermediate phase is separated by phase 
boundaries from the quantized phases, and the critical 
behavior at these transitions is governed by the nontriv- 
ial fixed points. At these fixed points, kxy is equal to 
the quantized value in the neighboring quantized phase. 
The critical exponents for these transitions are unknown 
at present. Experimentally, one would see plateaus in 
kxy, separated by intermediate regions, and the width 
of the latter will stay nonzero as the system size goes 
to infinity, and as the temperature goes to zero. In the 



intermediate regions, 



'^xy 



will vary continuously to in- 



terpolate its neighboring quantized values, and kxx will 
have a peak, the height of which will grow logarithmically 
with increasing system size or inverse temperature. We 
emphasize that the charge transport properties are still 
either superconducting or quantized Hall, depending on 
the system considered, and unaffected by the transition 
in the quasiparticles (there would be a collective mode 
(phonon-like) contribution to thermal transport in a neu- 
tral superfluid case, such as a He^ film). 

We should respond to one possible objection to our 
claim that there will generically be an intermediate 
metallic phase in the class D problem (this point is raised 
in Ref. |Q). This objection begins with the pure prob- 
lem, in which (since we assume noninteracting quasipar- 
ticles) the critical theory is a Majorana fermi field with 
a mass term that changes sign, and then considers weak 
disorder as a perturbation of this continuum theory. The 
similar pro blem of a Dirac field has been analyzed in re- 
cent years [|lOO| , and a central argument is that for E — 0, 
one can consider the problem as a Dirac field in two (Eu- 
clidean) dimensions. It is then argued that there are only 
a few possible random terms that can be included in the 
2D Dirac action that are marginal or relevant by power 
counting at the pure fixed point. These terms, which are 
bilinears in the Dirac field since the problem is noninter- 
acting, are a random U(l) vector potential and two types 
of random mass term. For the random Majorana fermion, 
no U(l) vector potential is possible, and there is a unique 
mass term. It is further known that the mass term is 
marginally irrelevant for weak disorder; this type of ran- 
domness arises when one considers the 2D Ising model 
with disorder in the intersite Ising couplings [101|. If one 
assumes that the disordered paired system we consider 
here must fall into this scheme, then the only possible 
randomness is irrelevant, and there will be a transition 
directly between the disordered versions of the weak- and 
strong-pairing phases, with the critical properties of the 



pure system, up to logarithmic corrections, in disagree- 
ment with our prediction of an intermediate phase and a 
different universality class. 

There is, however, a form of disorder not considered in 
this argument. In Sec. ^, we discussed vortices in the gap 
function, and in Subsec. [II Q related their description in 
terms of the gap function and the vector potential of the 
underlying problem to the vielbein and spin connection 
that appear in the general Majorana action. We should 
consider the possibility that these are random; indeed, 
the general analysis of AZ requires generic randomness, 
even though the Dirac or Majorana actions do not ex- 
plicitly appear there. Not all of the random fluctuations 
of these quantities are necessarily relevant. For exam- 
ple, small random fluctuations of the magnitude of the 
vielbein (i.e., of |A|) around its non-zero mean are ir- 
relevant by power counting. Also, if the gap function 
is nonzero the vector potential can mostly be gauged 
away because the superconductor is a Meissner phase. 
However, this is clearly not true at a vortex, and for 
general probability distributions of the disorder, vortices 
will be present. In the FQHE application, there is un- 
derlying potential disorder which couples to the particle 
density, and for unbounded distributions of disorder, or 
for sufficiently strong disorder with a bounded distribu- 
tion, it will nucleate vortices (FQHE quasiparticles) in 
the ground state, which can occur in isolation from other 
vortices since they have finite energy. It is not difficult 
to see, either intuitively or formally, that a small density 
of randomly placed vortices will be highly relevant at the 
massless Majorana critical point. Intuitively, they intro- 
duce destructive phase interference. Formally, in the 2D 
Majorana action, the effect is to insert the spin fields of 
the Majorana theory (so-called because they represent 
the Ising spins in the related Ising model), which (in the 
same gauge choice as earlier) cause square-root branch 
points in the Majorana field |l92| . Such random vortices 
do not seem to have been considered in previous work on 
random Dirac fields. They are relevant because, while 
the spin field in the critical Majorana theory has scaling 
dimension 1/8 (which corresponds to a relevant pertur- 
bation even in the non-random case), on replicating or 
supersymmetrizing the system, the spin fields act on all 
the components simultaneously, and hence their dimen- 
sion is then x 1/8 = 0. Thus the coefficient with which 
these fields appear in the action after averaging has scal- 
ing dimension 2, showing they are strongly relevant [P8[ . 
This will cause an RG flow away from the pure Majorana 
fixed point, and we expect the generic behavior of class 
D, with the intermediate metallic phase, to result. An- 
other possible form of disorder is that which violates the 
I — —1 symmetry, that is general p^ + ipy pairing with 
random coefficients. If both were completely random, the 
average would restore parity, that is the symmetry of re- 
fiection in any line, and prevent the existence of nonzero 
kxy, SO this form of disorder is not completely acceptable 
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in our physical systems; we must allow for a net violation 
of parity. Thus, the generic disorder that defines class D 
should include all of these relevant effects, and randomly 
placed, isolated vortices alone are a relevant perturbation 
that leads to a flow away from the pure Major ana criti- 
cal theory. Notice that, with random vortices included, 
the ordered phase we are describing is no longer a super- 
conductor in the strict sense, due to the random phases 
frozen into the gap function. Instead it is what has been 
called a gauge glass, since the order in the gap function 
is similar to that in a spin glass. 

In contrast, in the random bond Ising model (a subject 
also raised in Ref. ^^), negative bonds produce frus- 
trated plaquettes, and a string of negative bonds pro- 
duces two frustrated plaquettes at the ends of the string. 
The insertion of a semi-infinite string of negative bonds is 
the definition of the Ising disorder operator of Kadanoff 
and Ceva [102|, which is dual to, and at the critical point 
has the same scaling dimension as, the Ising spin field. 
However, in this case the analysis of weak randomness 
in terms of a random mass is justified, because for a low 
probability of negative bonds the disorder variables (vor- 
tices) appear only in closely-spaced pairs, not in isola- 
tion. In the continuum critical field theory, the operator 
product of two disorder operators, which represents the 
close pair, produces the Majorana mass operator, and 
thus the randomness generates the random mass term. 
We see that the distribution of the randomness in the 
Ising problem with weak bond disorder differs substan- 
tially from the problem we wish to consider. While the 
random bond Ising model in 2D does have a direct tran- 
sition between the two quantized k^y phases, and a mul- 
ticritical point, there is no reason to suppose that these 
occur in class D. 

The properties of the disordered version of the MR 
phase, which has a nonzero quantized kxy, are subtle. 
Since we have assumed that isolated vortices are possi- 
ble, which are localized FQHE quasiparticles analogous 
to those in the usual states on the plateaus in the FQHE, 
these will carry zero modes, and there will be 2"~^ many- 
particle states when the system contains 2n vortices. 
These are very nearly degenerate when the vortices are 
well-separated, since the energy splittings are expected to 
decrease exponentially in the separation of the vortices. 
Nonabelian statistics of additional vortex (FQHE quasi- 
particle) excitations should be considered in terms of ex- 
changes of such quasiparticles separated by many times 
the typical separation of the vortices in the ground state. 
Then the fcrmion zero modes of the ground state can 
interfere with those on the added quasiparticles, compli- 
cating the nonabelian statistics properties. Further study 
of these effects is beyond the scope of this paper. 

Finally, we return as promised to the case of the 
double-layer system with tunneling t. Then the phase 
boundaries at non-zero t between the weak-pairing, 
abelian phase, the MR phase, and the strong-pairing 



phase will be broadened and replaced by an intermedi- 
ate region in which metallic thermal properties will oc- 
cur, again with sharp phase boundaries between this and 
the other phases. As i ^ 0, this intermediate phase 
will shrink in width to become a single point at i = 0, 
where we have already explained that a direct transition 
in the IQHE universality class occurs. Thus this transi- 
tion broadens to become the intermediate metallic phase 
at finite t. We expect that at sufficiently small t, there is 
a single region of the metallic phase, which interpolates 
between k^y values differing by two steps, which are the 
values in the t = weak and strong-pairing phases. As t 
increases, a point is reached at which another plateau in 
Kxy appears, which is the MR-type phase. Such behavior 
is allowed by our flow diagram, if we plot it from 6 = to 
9 = 47r, and the initial values lie on an arc between the 
quantized fixed points at those 9 values that avoids the 
basin of attraction of the quantized fixed point a.t 9 ^ 2tt 
completely for small t, but not for larger t. 



VI. CONCLUSION 

In this paper, we have considered exotic properties of 
P and T-violating paired states of fermions in two di- 
mensions, and the relation to the FQHE using pairing of 
composite fermions. The results have been summarized 
in the Introduction. 

To conclude, our main results are: (i) p-wave pairing 
in spinless or spin-polarized fermions in the weak-pairing 
phase leads to the properties also found in the FQHE 
in the MR states, and supports the ideas of nonabelian 
statistics as a robust property, at least in the case of 
a pure system. Such statistics will also occur for the 
vortices in such a p-wave state in general in charged su- 
perfluids, and in neutral superfluids modulo U(l) phase 
factors that arise from the compressible charge sector; 
(ii) in an A-phase type p-wave phase, statistics may be 
abelian, though tunneling or a Zeeman term may lead 
to a transition to a MR phase; (iii) in the d-wave spin- 
singlet case, the HR state corresponds to the transition 
point, and, from now on, may be disregarded in consider- 
ing generic spin-singlet FQHE systems, which will most 
likely be in the weak-pairing phase. The latter is abelian, 
but has a nonzero Hall spin conductivity, and spin- 1/2 
chiral Dirac fermions on the edge |5^ ; (iv) disorder does 
not destroy the phases in question, but may modify the 
MR phase in an essential way. In the spinless p-wave 
case, randomly-placed vortices are a relevant perturba- 
tion of the pure transition, and there is an intermediate 
phase with metallic thermal conductivity properties due 
to delocalized BCS quasiparticles. 

Issues which remain to be addressed include the full 
effective theories for the states, and for the transitions 
between them, and the effect of interactions on the ran- 
dom problems. Also, a direct derivation of nonabelian 
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statistics in terms of pairing of fermions in the MR case 
is desirable. 

One further comment on tunneling into the edge is in 
order. Such tunneling could provide a useful diagnos- 
tic for the paired states in the FQHE. Since the fermion 
excitations on the edge in the weak-pairing phases are 
now always Dirac (or Majorana), their contribution to 
the exponent is always the same. Thus at filling factor 
V — \jq (where q will be 2, 4, . . . , for fermions such as 
electrons), the current will scale as / ~ V"" with a — q+\ 
in all the weak-pairing or MR phases. In contrast, in the 
Halperin-type paired states, we have a = iq [Q. The 
former result is t he sa me as in the compressible Fermi- 
liquid- like states 1 103 1 . Assuming that edge theories of 
fully-gapped bulk states must alwa ys b e unitary confor- 
mal field theories, as argued in Sec. IV, the exponent for 
tunneling into an edge on which all modes propagate in 
the same direction must in fact always be an odd integer, 
as a consequence of the Fermi statistics of the electrons. 
For the 5/2 state, the I ^ V contribution of the LLL will 
presumably dominate. 

Note Added: In view of a suggestion which has cir- 
culated, that the Fermi-liquid-like phase of Ref. ||2^ may 
generically have an instability to pairing in some channel, 
albeit at extremely low energy or temperature scales, we 
will consider here the case of a weak-pairing phase of spin- 
polarized fermions in an arbitrary angular momentum / 
eigenstate {I must be odd). Similar arguments as before 
show that there will be |^| chiral Majorana fermion modes 
on an edge, and correspondingly 2l'l"~^ degenerate states 
for 2n vortices. Since |^| must be odd, this always leads to 
nonabelian statistics, with the same monodromy proper- 
ties as for Z = —1 up to Berry phase factors (because each 
added pair of Majoranas makes a Dirac fermion which 
contributes only abelian effects). 



ACKNOWLEDGMENTS 

We thank M.P.A. Fisher, LA. Gruzberg, F.D.M. Hal- 
dane, B.I. Halperin, T.-L. Ho, A.W.W. Ludwig, M.V. 
Milovanovic, G. Moore, E.H. Rezayi, and T. Senthil for 
informative comments and helpful discussions. N.R. is 
grateful for the hospitality of the Institut Henri Poincare, 
Paris, where this paper was completed. This work was 
supported by NSF grants, numbers DMR-9157484 and 
DMR-9818259. 



APPENDIX A: QUANTIZED HALL 
CONDUCTIVITY FOR SPIN 

In this Appendix, we provide a detailed derivation of 
the Hall conductivity in (iso-) spin transport in the d- and 
p-wave (A-phase) cases. This is equivalent to showing 
that the induced action for the system in an external 



gauge field that couples to the (iso-) spin is a CS term. 
In the d-wave case, the system is spin-rotation invariant, 
so we obtain an SU(2) CS term, while in the p-wave 
case, there is only a U(l) symmetry, so we find a U(l) 
CS term. In both cases, the Hall spin conductivity is 
given in suitable units by a topological invariant. Within 
the BCS mean field approach, using a suitable conserving 
approximation, this topological invariant is the winding 
number of (wk, "yic) discussed in Sec. ||, and is therefore an 
integer, which is the statement of quantization. We argue 
that the quantization in terms of a topological invariant 
is more general than the approximation used. 

Considering first the spin-singlet paired states, we use 
the Nambu basis where the symmetries are transparent. 
Define 



* 



-( ' 



with 



(Al) 



(A2) 



so that $ transforms as a tensor product of particle-hole 
and spin-space spinors. We consider an interacting sys- 
tem and approximate it as in BCS theory, then with a 
minimal coupling to the gauge field, we use a conserving 
approximation to obtain the spin response. In Fourier 
space, we should note that 
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Ck 



(A3) 



In the Nambu basis, the kinetic energy becomes (again, 
K = H - fiN) 



^o^Yl ^klCkT^kT + Cki^ki) 
k 



(A4) 



where ^^ = \k\'^/{2m) — ^ is the kinetic energy, containing 
the bare mass tti, and the products in the spinor space 
are understood. Products like CTz / act on the Nambu 
spinors, with the first factor acting in the particle- hole 
factor, the second in the spin-space factor. The interac- 
tion term, for a spin-independent interaction V, is 



i^int, = 



I E ^1 ^ *k+q('^- ® /)*k*L_q(a. ® /)*k' 



kk'q 



(A5) 



Here the colons : . . . : denote normal ordering, that is 
all the c^s are brought to the left. In the BCS-extended 
Hartree-Fock approximation, the effective quasiparticle 
Hamiltonian (for later reference) is 
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KcH = XI *k Kk(fT2 ® /) + Re Ak(CT:. ® /) 

k 

-ImAk(cry«)/)]*k- 



(A6) 



This is for singlet pairing, where A„k = Ak, and not just 
for d-wave. Here fk is Ck P^^^ ^^^ Hartree-Fock correc- 
tions. If we define a vector 

Ek = (ReAk,-IniAk,a) (A7) 

then the quasiparticle energy Ey^ = |Ek|j and 

Xeff-X<(Ek-^®/)*k. (A8) 

k 

In the Nambu notation, it is clear that K — Ko + Kint, 
and KcS, are invariant under global SU(2) rotations that 
act on the spin-space, that is the second factor in the 
tensor products. The spin density, the integral of which 
over all space is the total spin and generates such global 
transformations, and the spin current densities are given 
by 

^0 (q) = ^ E *k-q/2(^ ® ^a)*k+q/2 (A9) 

k 

^r(q) = ^ E ^<-q/2('^- ® ^a)*kW2, (AlO) 

k 

where i = x, y is the spatial index, and a — x, y, z 
is the spin-space index. Spin conservation implies the 
continuity equation, as an operator equation, 



dJVdx^ = 0, 



(All) 



where /i = 0, a;, y, and the summation convention is in 
force. 

So far we have not introduced a gauge field for spin. 
Since the spin is conserved locally, we can turn the sym- 
metry into a local gauge symmetry by introducing an 
SU(2) vector potential, and making all derivatives co- 
variant. The effect on K is to add the integral of 
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A" J" + ^AlAl^ii^a, I)^. 
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(A12) 



The gauge field is to be used solely as an external source, 
with which to probe the spin response of the system, and 
then set to zero. 

If we now consider integrating out the fermions, then 
we can obtain an action in the external gauge fields, 
which can be expanded in powers of A" The zeroth- 
order term is the free energy density, times the volume 
of spacetime, and the first-order term vanishes by spin- 
rotation invariance. The second-order term corresponds 
to linear response: the second functional derivative with 
respect to A" at A"^ = 0, yields the (matrix of) linear 
response functions. In particular, the spatial components 



yield the conductivity tensor in the usual way. Therefore 
we consider the imaginary-time time-ordered function. 



n 



ab 



<J:.{q)4{-q)), 



(A13) 



where time-ordering is understood, and from here on we 
use a convention that p, q, etc stand for three-vectors 
p = {po , p) , and further po — iuj is imaginary for imagi- 
nary time. For fi = v = i — x or y, an additional "dia- 
magnetic" term nS"''' /Am is present in 11°^, which we do 
not show explicitly. As consequences of the continuity 
equation and the related gauge invariance, 11^;^ must be 
divergenceless on both variables, g^II^i/ = qi^Hni, = 0. 
To maintain these when using the BCS-Hartree-Fock ap- 
proximation for the equilibrium properties, one must use 
a conserving approximation for the response function, 
which in this case means summing ladder diagrams (com- 
pare the charge case in Ref. B, pp. 224-237). 

One begins with the BCS-Hartree-Fock approximation, 
which can be written in terms of Green's functions as (we 
consider only zero temperature, and J dpo is along the 
imaginary po axis throughout) 



G-\p)=po 
E(p) = 



_fk_ 

(27r)3 



®/-S(p), 



(A14) 
'I)V{k-q). (A15) 



Note that G{p) and S(p) are matrices acting on the ten- 
sor product space. The equations are solved by 



G \p)^po-'Ep-(T®I, 



(A16) 



(we write 1 for / (g) /) as one can also see from the effec- 
tive quasiparticle Hamiltonian Kcs, and Ap obeys the 
standard gap equation. 

In the response function, the ladder series can be 
summed and included by dressing one vertex, to obtain 
(again not showing the diamagnetic term) 



KM 



d?p 



tr[7;(p,p + <7)G(p + g) 



xTl{p + q,p)G{p)\, 



where 7"^ is the bare vertex. 



Zm 
and F" is the dressed vertex satisfying 



(A17) 

(A18) 
(A19) 



Tl{p + q,p)=-,lip + q,p)+i 



(fk 

(2^ 



a:, (g) IG{k + q) 



xTl{k + q,k)G{k)a^(g)IV{p-k). (A20) 

At small q, we can obtain useful information about this 
function from the Ward identity that results from the 
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continuity equation. The particular Ward identity we 
use here is an exact relation of the vertex function to 
the self-energy, and the conserving approximation (the 
ladder series) was constructed to ensure that it holds also 
for the approximated vertex and self energy functions. 

Following Schrieffer's treatment M, we consider the 
vertex function with external legs included: 



A^(ri,r2,r3) = (J„n^3)*(ri)*t(r2)), 



(A21) 



for spacetime coordinates ri, r2, r^. Applying d/dr^f^ to 
both sides and using the operator continuity equation, 
we obtain the exact identity in Fourier space 

q^^V^ip + q,?) - -I®aaG-\p) - -G-\p + q)I®aa. 

(A22) 

Since G^^ is trivial in the spin-space indices, it commutes 
with / (g) (Tq. Hence at q ^ 0, the right-hand side van- 
ishes, so T{p -\- q,p) has no singularities as q ^ 0. This 
differs from the charge case, for example, where this cal- 
culation (using the ladder series approximation) leads to 
the discovery of the collective mode p^ . Since the spin 
symmetry is unbroken, no collective mode is necessary to 
restore this conservation law, and so there is no singular- 
ity in the vertex function for spin. 

One can verify that the Ward identity is satisfied using 
the BCS-Hartree-Fock G"^ and the ladder series for F. 
At q = 0, this yields the important results 



1 



T-^Xp,p) = -,^d^G'\p)I®aa. 
or explicitly. 



(A23) 



^o{p,p)^l^I®<ya, 

T1 ^ ^]^d,G-\p)I ® aa, 



(A24) 



where di and 9^ stand for d/dpi^ d jdp^^ from here on, and 
the extra minus in the first relation is consistent because 
implicitly g^F^ ^qa^^- qi^i- 

We now calculate 11 at small q. To zeroth order, use 
of the Ward identity shows that the J-J function gives 
zero, except when [i = v ^ i. In that case, it reduces to a 
constant that cancels the diamagnetic term also present 
in just that case. Hence we require only the part first- 
order in q. In the expression for H above, we first shift 
p — > p — iq, so that q no longer appears in any bare ver- 
tices, but does appear in the Green's functions on both 
sides of the ladder, between the rungs which are the in- 
teraction lines. Hence to first order, we obtain a factor 
±ii9^G = '^\GdG~^G in place of G in one position in 
the ladder. Since there may be any number of rungs 
(including zero) between this and either of the vertices 



at the ends, the terms can be summed up into a lad- 
der dressing each vertex, evaluated at g = 0. Hence we 
obtain to first order 



KM 



Sp 



lx^1Xp,p)q^d^GYl{p,p)G{p) 

(A25) 



2 7 (27r)3 L ^i 

-Vl{p,p)G{p)Vl{p,p)qy.dxG\ 

Using the Ward identity, this becomes 
Sp 



KM 



-,q\ 



ic, ..tr {{I ® aa){I ® (7b)Gdf,G ^ 



[GdxG-\Gd^G-^]} 



(A26) 



Since the G's are independent of the spin-space indices, 
the explicit cr's factor off, and the result is (5°^ times a 
spin-independent part. The latter can be simplified using 
the BCS-Hartree-Fock form of G, by writing the latter 
as 



Gip) 



po + Ep ■ (T 1^ I 
Vl-El ■ 



(A27) 



The spin-independent factor contains e^j/A since it is an- 
tisymmetric in these labels. Keeping track of the signs, 
we find for the quadratic term in the induced action 



with M. given by the topological invariant 



(A28) 



(A29) 



The right hand side is exactly the winding number ni 
discussed in Sec. ||, and is an integer as long as E is a 
continuous, differentiable function of p; it is 2 for the 
d-wave case. 

To ensure SU(2) gauge invariance, the CS term should 
include also a term cubic in A, with no derivatives. For 
this term we evaluate the triangle one-loop diagrams with 
three insertions of J, with each vertex dressed by the 
ladder series. Setting the external momenta to zero, the 
Ward identity can be used for all three vertices, and the 
result can be seen to be 

n^tA(0, 0) ^ - ^ / (03 *>■ [(^ ® <^a)Gd,G-^ 

X {{I ® ab)Gd,G-\iI ® a,)GdxG-^}] . 

(A30) 

The anticommutator { , } arises since the result must be 
symmetric under permutations of the index pairs /x, a, 
etc. The product cra<Tb<^c: when traced over the spin- 
space indices, yields a factor 2ieabci which is antisymmet- 
ric, and so the remainder must contain e^^A to maintain 
symmetry; the rest of the structure is the same as be- 
fore. Hence the full result is the SU(2) CS term, which 
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we write in terms of the 2x2 matrix vector potentials 
in- A°- 



A., 



An 



(fxef,^xtT{A^d^Ax + -Af^Ar^Ax). 



(A31) 



Here k is the conventional notation for the coefficient of 
such a term, in this same normalization; if we wished to 
quantize the theory by functionally integrating over A, 
we would need fc = an integer. In our case k = Ai/2 — 1 
for d-wave. 

For the spin-triplet case with an unbroken U(l) sym- 
metry, we must use the fact that A_k = — Ak- For exam- 
ple, in the two-dimensional A-phase, as occurs in the 331 
state in the double-layer FQHE system with zero tun- 



neling, discussed in Sec. Ill, the pairs are in the isospin 
Sz = triplet state tiij + Ulj, and the U(l) symmetry 
generated by Sz is unbroken; we recall that the under- 
lying Hamiltonian is not assumed to have a full SU(2) 
symmetry. The effective quasiparticle Hamiltonian ( |3§| ) 
becomes, in the Nambu-style notation, 

K^S = 51 *k Kk(fT^ ® /) + Re Ak(CT:. ® CTz) 



ImAkK0cr^)]^k- 



(A32) 



The U(l) vector potential A^ couples to Sz, and the ver- 
tex functions contain I ® Gz, which commutes with the 
BCS-Hartree-Fock Green's function G. The tensors ap- 
pearing in the three terms in K^s obey the same algebra 
as the three in that for the spin-singlet case (where they 
were trivial in the second factor), and as in that case 
commute with I ® Uz- Consequently, the derivation for 
the induced action to quadratic order in A^ is similar 
to that for the SU(2) singlet case above, and the traces 
in the Nambu indices can be carried out with the same 
result as before, to obtain the abelian CS term 



An 



d^rA, 



dA, 
drx 



^fivX 1 



(A33) 



and no cubic term. In this case, M is again given by the 
winding number m which is or ±1 in the p-wave strong 
and weak-pairing phases (respectively) discussed in this 
paper. 

We note that the effect of the vertex corrections we in- 
cluded as ladder series is to renormalize the g = vertices 
as shown in eq. ( A23| ) for the spin-singlet case, and use 
these in one-loop diagrams with no further corrections. 
This corresponds to the minimal coupling p — > p — A 
in the action, as one would expect by gauge invariance. 
If we assume such a coupling, and treat the low-energy, 
long-wavelength theory near the weak-strong transition 
as Dirac fermions with rclativistic dispersion and min- 
imal coupling to the external gauge field, then the ex- 
pression for Ai as an integral over p covers only half the 
sphere in n space, and we would get ±1 (d-wave), ±1/2 



(p-wave). The missing part results from the ultraviolet 
regulator in the field theory version of the calculation 
p2[ , or from a second fermion with a fixed mass in some 
lattice models [100|. In our calculation, the remainder 
is provided by the ultraviolet region, where Ak — *■ as 
k ^ CX3. At the transition, /i = 0, the map is discontinu- 
ous and covers exactly half the sphere in the p-wave case, 
so Al = 1/2, as in other problems. In the d-wave case 
with rotational symmetry, the value of |uk/uk| as k — > 
is nonuniversal, as noted in Sec. IVA, and hence so is 



the value of a^. at the transition. This is a consequence 
of the non-relativistic form of the dispersion relation of 
the low-energy fermions in this case. We may also note 
that for a paired system on a lattice, as in models of 
high Tc superconductors, a similar calculation will give 
an integral over the Brillouin zone, which is a torus T^, 
instead of the k plane which can be compactified to S^ . 
But maps from T^ to S^ are again classified by the in- 
tegers, and the integer winding number is given by the 
same expression, so quantization is unaffected. 

We can also argue that the quantization result away 
from a transition is exact in a translationally-invariant 
system, at least in all orders in perturbation theory. For 
this we use the form in eq. (A26) or ( A3C| ), where the 
Ward identity for the vertex has been used. Diagram- 
matically, it is clear that the exact expression can be 
similarly written, using the exact (i.e., all orders in per- 
turbation) Green's function and vertex function. (This is 
also true when the CS gauge field interaction is included.) 
The Ward identity that relates them is exact, and the re- 
sult for cr| is of the same form as shown. The next step, 
the frequency integrals, cannot be done explicitly in this 
case, because the precise form of the Green's function is 
unknown, and the analogs of ^k, Ak (or of Wk, fk) do not 
exist. The latter do not exist because in general the poles 
in the Green's function, which would represent the quasi- 
particles, are broadened by scattering processes, except 
for the lowest en ergie s for kinematical reasons. However, 
the form in eq. ( A26 ) is itself a topological invariant, as 
we will now argue. As long as there is a gap in the sup- 
port of the spectral function of G, G{p) is continuous and 
differentiable on the imaginary frequency axis, and tends 
to I®I Ipq as po ^ ±ioo. Thus G~^ exists and never van- 
ishes. Considering the spin-singlet case for convenience, 
the spin-space structure is trivial, so we may perform the 
corresponding traces, and then G or G^^ is a 2 x 2 ma- 
trix, with the same reality properties on the imaginary 
Po axis as in the BGS-Hartree-Fock approximation. (The 
spin-triplet case should work out similarly, because of 
the algebraic structure already mentioned.) It thus rep- 
resents a real non-zero 4-component vector, in R^ — 0, 
which topologically is the same as S^ . S^ is obtained by 
dividing G by its norm, (trG^G)^'^, and the normalized 
G is a 2 X 2 unitary matrix with determinant —1, so it 
lies in S^ . The k space can be compactified to S'^ as 
before, and the frequency variable can be viewed as an 
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element of the interval I = (—1, 1), so the integral is over 
S'^ X I. However, since the limit of the Green's function 
as po — > zLioo for fixed k is independent of k, we can 
view this as simply S^. Thus we are dealing with maps 
from S^ to S^ , the equivalence classes of which are clas- 
sified by the homotopy group 7r3(S'^) — Z. The integral 
we have obtained simply calculates the integer winding 
number or Pontryagin index of the map, when properly 
normalized {G can be normalized to lie in SU(2) without 
affecting the integral). This establishes the quantization 
of a^ in a translationally-invariant system with a gap, 
at least to all orders in perturbation theory, and prob- 
ably can be made fully non-perturbative (as the Ward 
identity is already). 
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